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. BACKGROUND  AND  SUMMARY  OF  FINDINGS 


Background  | 

Venting  systems  are  employed  on  marine  vessels  to  relieve  pressure  dif- 
ferences which  arise  between  the  cargo  tank  and  the  ambient.  Based  on  cargo 
flammability  and  toxicity,  either  of  two  types  of  venting  systems  is  selec- 
ted: 

• Low  Flammability  and  Toxicity 

Individual  gooseneck  vents  fitted  with  flfime  screens  and  ball- 
check  or  pressure-relief  valves:  Here  the  pressure  build-up  dur- 
ing loading  is  negligible  because  open  gauging  is  permitted. 

• High  Flammability  or  Toxicity  j 

Tanks  manifolded  to  vent  risers  fitted  with  PV  valves  and/or  flame  | 

: arrestors:  Here  the  normal  pressure  build-up  during  loading  is  < 

0-3  psig  (depending  on  PV  settings)  because  restricted  or  closed  j 

gauging  is  required.  I ^ 

1 J 

In  the  latter  case  if  the  cargo  pumping  rate  exceeds  the  vapor  relief  capa- 
city, or  if  there  is  an  accidental  overfill,  the  pressure  difference  can 
rise  to  several  atmospheres.  Depending  upon  the  design  of  the  bulkhead 
structure,  tank  failure  may  occur  in  such  "overpressure"  cases.  The 
overpressure  hazard  and  the  criteria  for  controlling  it  are  being  analyzed 
by  Arthur  D.  Little,  Inc.  (ADL)  as  part  of  contract  No.  DOT-CG-42,357-A. 

Based  on  this  study,  we  have  formulated  a procedure  to  evaluate  the  over- 
pressure of  any  given  cargo  transfer  operation,  based  on  cargo  properties, 
loading  rate,  tank  characteristics,  and  vent  system.  Design  and  opera- 
tional guidelines  are  suggested  to  avoid  overpressurization.  The  procedure 
and  guidelines  are  described  in  this  report,  along  with  the  engineering 
' analyses  which  support  them. 

Summary  of  Findings 

Based  on  a review  of  ABS  rules  and  structural  analyses  of  decks  and 
bulkheads  of  specific  barges  and  tankers,  we  estimate  that  failure  of  cargo 
tanks  occurs  at  pressures  exceeding  about  6 psig. 

Typical  venting  systems  in  use  today  have  more  than  adequate  capacity 
for  gas  venting  but  Inadequate  capacity  for  liquid  overfill.  For  example, 
f the  predicted  pressure  rise  during  gas  venting  does  not  exceed  6 psig 
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unless  (a)  the  effective  L/D  of  the  vent  system  Is  about  6,000  (more 
than  ten  times  the  normal  value)  or  (b)  the  loading  rate  Is  three  times 
normal.  The  combinations  of  loading  rate  and  effective  L/D  which  lead 
to  hazardous  tank  pressures  are  given  In  Table  IV-2  (page  41) . Cargo 
volatility  has  a relatively  small  effect  on  the  hazard  evaluation;  a 
typical  volatile  cargo  adding  less  than  2 pslg  to  the  normal  pressure 
differential. 

An  analysis  of  the  pressure  buildup  during  overfill  suggests  that 

O 

the  maximum  permissible  ratio  of  loading  rate  (ft  /sec)  to  vent  cross- 
2 

sectional  area  (ft  ) Is  9 ft/sec.  Many  cargo  transfer  operations 
currently  equal  or  exceed  this  value,  having  ratios  of  30  ft/sec  and 
reater  In  order  to  meet  reasonable  vessel  turnaround  times.  Further- 
e,  tank  failure  Is  expected  within  a very  short  time  (less  than 
iO-seconds)  after  the  tank  becomes  liquid  full.  Therefore,  It  appears 
that  the  liquid  overfill  situation  must  be  prevented  before  It  occurs 
rather  than  taking  action  after  It  occurs. 
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A.  Tank  Configurations 

The  cargo  space  of  chemical  tankers  and  barges  Is  partitioned  into 
10  to  30  individual  tanks,  both  for  structural  reasons  and  in  order  to 
segregate  different  cargoes.  Ranging  from  500  to  70,000  DWT  in  total  dis- 
placement and  with  a typical  size  in  the  range  35,000  to  40,000  DWT,  "drug- 
store" tankers  are  an  order  of  magnitude  lower  in  capacity  than  crude  oil 
carriers.  A typical  individual  tankshlp  tank  volume  is  about  12,000  bbl 
or  60,000  ft^.  About  50  chemical  tankers  of  under  10,000  DWT  are  currently 
being  built  (Farrell  [1974]).  The  tank  and  piping  layout  of  a typical  chemi 
cal  tanker  is  shown' in  Figure  II-l. 


Figure  Il-l  - Typical  Cargo  Segregation  Arrangement 
for  Chemical  Tanker 


(Farrell  [1974]) 


In  addition  to  the  200  billion  ton-miles  of  yearly  U.  S.  freight  trans- 
port by  tanker,  an  equal  or  greater  amount  of  transport  was  done  by  about 
2,500  barges  used  on  Great  Lakes  and  Inland  waterways,  with  a capacity  rang- 
ing from  7,000  to  40,000  bbls.  An  Individual  tank  on  a barge  would  nominally 
be  about  2,000  bbls  or  11,000  ft^.  The  tank  layout  of  a typical  250  x 50  x 12 
ft  barge  Is  shown  In  Figure  Ti-2. 
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Each  Tank  2,700  bbls  or  15,000  Cu.  Ft. 

Barge  Dimensions  250  x 50  x 1 2 Ft. 

FIGURE  11-2  TANK  CONFIGURATION  OF  BARGE 


B.  Cargo  Transfer  Procedures 

Normally,  the  tanks,  piping,  and  pumps  are  matched  and  sized  to  com- 
plete loading  or  offloading  of  the  entire  vessel  within  approximately  5 
to  10  hours.  The  tanks  are  loaded  by  subgroups  of  about  four  tanks  to  each 
piping  manifold.  Based  on  the  tank  sizes  noted  In  Section  II-A,  the  cargo 
transfer  rates  are  normally  about  1,500  bbl/hr  per  tank  for  a 12,000  bbl  tank 
on  a tankshlp  and  2,500  bbl/hr  for  a 20,000  bbl  barge.  The  pumping  rate  for 
a given  vessel  size  Is  shown  In  Figure  II-3. 

For  most  tankers  and  barges,  the  loading  pumps  are  on-shore,  but  off- 
loading pumps  are  part  of  the  vessel's  equipment.  In  both  cases,  the  cargo 
transfer  operation  requires  adequate  communication  between  vessel  crewmen 
and  shore  pump  operators.  Recently,  a move  toward  automation  of  dockside 
operations  has  been  clear,  with  Increased  responsibility  on  the  vessel 
crew.  In  some  operations  the  cargo  supply  hose  Is  connected  sequentially 
to  two  or  more  cargo  transfer  headers. 
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100  200 
Deadweight  (OOO's  of  tons) 


FIGURE  11-3  PUMPING  RATES  FOR  TANKERS 
(Paterson  and  Watson  [1974] ) 


c. 

Key  Cargo  Properties 

The  cargoes  carried  in  bulk  can 

be  divided  into 

six  major  classes 

(Farrell  [1974]): 

• 

Paraffin-based 

(Aliphatic) 

ethylene 
vinyl  chloride 
formaldehyde 
propylene 

methane 

acrylonitrile 

butadiene 

naphtha 

acetaldehyde 

• 

1 

Aromatic-based  i 

( 

t coal  tar 
naphthalene 
benzene 

toluene 

xylene 

styrene 

• 

Alcohol-based 

1 

\ ethanol 
j methanol 

propanol 
butyl  alcohol 

• 

Ester-based 

( ethyl  acetate 
< vfnyl  acetate 
' butyJ  acetate 

• 

Inorganics 

1 

f acids 

1 soda 

1 sulphur 

• 

Mixtures 

1 crude  oil 
( gasoline 

1 mineral  spirits 

kerosine 
lubricating  oils 

Two  properties  relevant  to  the  hazards  of  excessive  vapor  displace- 
ment and  liquid  overfill  are  the  vapor  pressure  and  viscosity,  respectively. 
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In  Table  II-l  are  collected  the  vapor  pressure  and  liquid  viscosity  of 
many  bulk  cargoes  of  concern. 
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Table  II- 1 


Properties  of  Bulk  Cargoes  1 


Vapor 

Liquid 

Pressure 

Viscosity 

Chemical 

(Psla  @ 100  F) 

(cp  @ 100  F) 

Acetonitrile 

3 

- 

Allyl  Chloride 

11 

- 

Benzene 

3.2 

.5 

t-Butyl  Alcohol 

1.8 

1.8 

i-Butyraldehyde 

6 

.35 

Cyclohexane 

3 

.7 

Diethylamine 

7.7 

,3 

Diisobutylene 

1.7 

- 

Distillate  (Straight 

Run)  1. 7 

.35 

Ethyl  Acetate 

3.1 

.37 

Ethyl  Alcohol 

2.1 

- 

Ethyl  Benzene 

0.4 

.54 

Ethyl  Ether 

18 

.2 

Ethylene  Glycol  Dimethyl  Ether  2.2 

- 

Gasoline 

Boils  at  100  to  400  F 

.35 

Heptane 

1.7 

- 

Hexane 

5 

.30 

Hexene 

6 

.22 

Isoamyl  Alcohol 

0.2 

2.6 

Isopropyl  Alcohol 

1.8 

- 

Isohexane 

7 

.25 

Isopentane 

Boils  at  82.2  F 

.23  @ 70  F 

Isoprene 

Boils  at  93.4  F 

.2  @ 80  F 

Isopropyl  Acetate 

2 

.41 

Methly  Acrylate 

3 

.4 

Methyl  Alcohol 

4.5 

Methyl  Ethyl  Ketone 

3 

- 

Methyl  Methacrylate 

1.5 

.45 

Mineral  Spirits 

0.1 

2.6 

Naphtha 

0.1 

2.6 

Oil-Penetrating 

0.1 

2.8 

Pentane 

Bolls  at  98.0  F 

.23  0 78  F 

Pentene 

Boils  at  85.8  F 

.20  0 68  F 

Propionaldehyde 

6.7 

.30 

Toluene 

1 

.48 

Vinyl  Acetate 

3.7 

.37 

■IMM 


D.  Vent  System  and  Constrictions 
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Two  types  of  vent  systems  are  found  on  cargo  vessels;  these  are  termed 
masthead  and  standpipe  systems.  In  the  masthead  vent  system,  shown  in 
Figure  II-4,  groups  of  tanks  are  manifolded  to  a common  header  which  runs 
along  the  deck  and  up  the  kingpost  or  mast  to  a height  which  may  exceed  1/3 
the  vessel  beam.  Spill  valves  may  be  Installed  where  the  smaller  pipes 
join  the  header  and  at  the  base  of  the  vertical  riser  to  prevent  overfilling 
to  the  top  of  the  outlet.  A flame  arrestor  or  PV  valve  may  be  placed  in 
the  vertical  riser.  The  masthead  system  is  more  prevalent  for  older  vessels, 
more  effective  for  dispersing  cargo  vapors,  less  economic  for  larger  vessels, 
and  is  required  by  Regulation  46  CFR  32.55-20  (b)  for  Grade  A liquid  cargoes. 

The  standpipe  system,  illustrated  in  Figure  II-5,  consists  of  a verti- 
cal vent  pipe  above  each  tank  designed  to  release  excessive  vapor  in  mix- 
ture at  a height  of  3 to  8 ft  above  deck  level.  The  pipe  is  normally  fitted 
with  a PV  relief  valve,  a flame  control  device,  and  a gooseneck  to  protect 
against  weather.  The  standpipe  system  is  common  on  newer  tankers,  larger 
tankers,  and  almost  all  barges. 

Both  vent  systems  present  a finite  constriction  to  the  gas  efflux  when 
the  P-V  valve  is  open.  The  vent  pipe  diameters,  which  range  from  2-1/2"  to 
12",  are  sized  according  to  the  anticipated  displacement  rates,  which  in 
turn  are  limited  by  the  size  of  the  cargo  filling  lines  due  to  the  resis- 
tance of  liquid  inflow.  The  resistance  to  gas  outflow  is  several  orders  of 
magnitude  less  than  the  liquid  pressure  drop  for  a given  vent  system.  Cer- 
tain vent  piping  systems  appear  to  have  been  sized  large  enough  to  accomodate 
a liquid  cargo  overfill;  these  systems  are  "overdesigned"  for  gas  venting. 
However,  as  we  shall  point  out  below,  for  most  vent  systems  the  capacity 
of  the  P-V  valves  and  other  vent  fittings  generally  are  marginal  for  li- 
quid flow,  giving  rise  to  pressure  drops  of  several  atmospheres  during  an 
overfill.  That  is,  many  vent  systems  appear  to  have  excess  capacity  for 
gas  venting  and  Inadequate  capacity  for  cargo  overfill. 

The  effective  length-to-dlameter  ratio,  L/D,  is  a convenient  index  of 
the  constriction  presented  by  a vent  system  to  either  liquid  or  vapor  flow. 

In  pressure  drop  calculations,  the  L/D  is  multiplied  by  the  friction  factor 
f;  the  latter  accounting  for  variations  in  flow  rate  and  wall  roughness. 

We  have  tabulated  in  Table  II-2  the  appropriate  L/D  for  various  pipe  fittings 
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Table  II-2 


Representative  Equivalent  Length  in  Pipe  Diameters  (L/D) 
of  Various  Valves  and  Fittings 


Standard  Valves 

L/D 

Globe  valves  (fully  open) 

340 

Gate  valves  (fully  open) 

13 

Gate  valves  (1/2  open) 

160 

Swing  check  valves  (fully  open) 

135 

Check  valves  (fully  open) 

150 

Butterfly  valves  (fully  open) 

20 

Standard  Fittings 

90  degree  standard  elbow 

25 

45  degree  standard  elbow 

16 

180  degree  "U"  or  gooseneck 

40 

Standard  tee 

Flow  through  run 

20 

Flow  through  branch 

60 

Vent  System  Components 

Flame  arrestor  (typical) 

200  + 100 

P-V  relief  valve  (typical) 

100  + 50 

Sources:  Data  on  standard  fittings  taken  from  Crane  Company,  Technical 
Paper  410. 

Data  on  USCG-approved  valves 
GPE,  Protectoseal,  and  Varet 

and  flame  arrestors  taken  from 
catalogues. 

and  valves.  For  a masthead  system,  the  effective  L/D  may  reach  several  hun- 
dred. This  is  indicated  in  Figure  II-6,  which  is  based  on  the  vent  system 
of  a typical  70,000  DWT  tanker. 


lote;  This  figure  is  for  illustrative 

purposes  only  and  is  not  intended 
to  depict  an  actual  arrangement. 


L/D  = 60  for  Tee 


Ullage  Opening 


1 2"  Dia  Header 


L/D  = 40  for  Gooseneck 


2 psi  Overfill  Relief 


Flame  Arrestor 
L/D  = 200 


L/D  = 25  for  elbow 


Total  Effactiva  Length  to  Diamatar  Ratio: 


6"  Dia 


P-V  Valve  (open) 
L/D*  100 


L/D  = 25  for  Elbow 


Item 

(itatniied  from  tank  outlet) 

Elbow 
Run  10.5  ft 
Elbow 
P-V  Value 


L/D -171  +473-644 

Effactiva  LfD  of  Effactiva  L/D 
Branch  Lina  of  Header 


E.  Normal  and  Extreme  Cargo  Transfer  Situations 


In  the  above  sections  we  have  described  the  tank  configuration,  the 
cargo  transfer  procedures,  the  cargoes,  and  the  vent  systems.  Under  normal 
situations  the..  Is  no  overfill,  and  the  vent  releases  vapor/alr  mixture 
at  100  to  130Z  he  volumetric  fill  rate,  with  a tank  pressure  rise  In 
the  range  of  0 to  pel  above  ambient.  The  normal  efflux  velocity  under 
these  conditions  Is  In  the  range  of  10  to  40  ft/sec  (Paterson  [1974]). 

Overpressure  can  occur  either  by  (a)  Inadvertently  overfilling  the 
tank  and  attempting  to  force  liquid  cargo  through  the  vent  system  or  (b) 
using  fill  rates  which  are  excessive  for  the  finite  constriction  of  a given 
vent  system.  The  parameter  values  of  these  two  extreme  cargo  transfer 
situations  are  contrasted  to  the  normal  values  below  In  Table  II-3. 

Overfill  Is  a widely  recognized  problem  and  was  cited  as  a contri- 
buting factor  In  the  gasoline  fire  and  explosions  of  tank  barge  "Ocean 
80"  In  Carteret,  New  Jersey,  In  1972  (USCG/NTSB-MAR-75-3) . In  this  case, 
however,  overfill  did  not  cause  overpressurization.  Apparently  this  par- 
ticular barge  was  designed  to  accomodate  overfill.  The  tolerance  for  error 
on  the  part  of  the  tankermen  or  dockside  pump  operators  Is  very  narrow  In- 
deed, without  improved  gauging  devices,  high-level  alarms,  spill  valves,  or 
overfill  tanks. 


Parameter  Values  for  Normal  and  Extreme  Transfer  of  Can 


o *o  0) 


lU 

M 

O 

CO 

O 

o 

> 

4J 

o 

•H  <U 

M-l 

* 

(0  U 

m 

in 

(0  a 

0)  (d 

1 

O 

o 

1 

U K 
X 

rH 

CM 

o 

H 

1 

o 

40 

m 

o 

o 

fo  I ^ a 

K cd  3 

I W 

4J 

O O 
m • <M  }-i 
iH  Q) 

0 0)  J3 
Cl.  3 4J 
Cl  tio  -a  cd 
J3  M 

(1) 

CO  O .O  60 


4J  in 

c 

u 

0 

iw 

4J  *H 

J2  >. 

iH 

-d  *0 

^ 4J 

ca 

O 

60  (d 

cn 

O 1 

•H  0 

a 

hj 

0 -H 

0 

•'O 

iH 

c 

CM  m 

r-l  ^ 

0 IV 

(0 

0 ^ 

p 

0) 

iH 

P 

4J 

0 

<0 

CM*^ 

fH 

CO  44 

Q 

0) 

0) 

•H 

•cn 

0) 

00 

H 44 

M 

> 

00 

(0 

Pd 

CO 

iH  4J 

> 

P 

rH  0) 

P 

0 

0 

Pd 

44 

0 

•H 

Pd 

• 

CO 

•H 

4J 

0 

0 

U 

Xi 

CO 

0 

P 

p 

'P 

M-l 

in 

r-l 

0 

0 

C«  0 

(d 

0)  0) 

•H 

bO 

Vi 

0 

f-H 

0 

U 

4J  CO 

0 

0 

1 

1 

iH 

-p  0 

0) 

0 

VO 

14 

0 

O .-H 
d M 
O CD 


M 

0 

0 > 

rH 

•H 

p.  0 

P > 

Vi 

Q> 

•H 

fH 

Q 

PdO 

44 

0 

Pd  0 

Vi 

44 

VI 

44 

M 

0)  P 

a 0 

3 13 

a VI 

3 dg 

Vi 

Vi 

> ^ 

p 

p 0 

P 

CC  iH 

p p 

bo  P 

00 

0 

S 

P 

P 

Vi 

X ^ 

> 

0 

X •H 

•H  p 

c 

Pd 

3 

^ Vi 

44 

P 

p 0 

•H 

60>CI 

p 

P4  43 

•H 

P 

Vi 

3 

^ . 

P 

VO 

P 

3 "H 

p 

VJ 

c 

P 

P 

rH 

•H 

P 

Vi 

P 

0 e 

V4  VI 

P P 

c 

P 

P4 

0 

Pd 

Vi 

P 

0 Vi 

P 

M a 

0 p 

P P 

p 

Vi 

> 

H 

0 

VI 

Vi 

Ev  0 

U 

9 0 

04  Vi 

CO  a 

Fv 

H 

a 

VI  44 

Q 

VJ 

•H 

Pd 

•H 

C 

0 

P 

Vi 

a 

0)  0 

I-) 

p 

0 

Vi 

3 

> 

44 

0 

0 

P 

X 

0 

•o 

Vi 

p 

n 

Pd 

Ni^ 

iH 

0 

TD  ti 

> 

p 

•H 

p 

•H 

> 

•H  P 

•H 

0 

> 

> 

44 

a a 

Vi 

•H 

• • 

M 

P s 

0 

Vi 

0 

0 

P 

M 

P 

•H 

P 

00 

00 

P 

c 

P 

> 2 

44 

•H 

VI 

VI 

VI 

P 

VI 

p > 

44 

Vi 

p 

p 

0 

H 

M 

Ex3 

u 

u 

iH  tH 

44 

VI  * 

VI 

0 

• 

•H 

P P 

P 00 

0 

m 

44 

(*4 

44  1 

44 

P U 

P rH 

VI  r-- 

0 

44 

0 

0 

P 

0 

P -3* 

0 

CM 

0 

^ 0 

P 0 

P P 

P 

• 

00 

00 

PO 

u 

P 

0 P 

> VI 

rH 

PO 

0 

> VI 

P > 

P 

it 

T3 

VI 

VI 

iH  P 

fH  dH 

P 0 

• 

VO 

0 

P 

Pd  00 

a X 

4-1 

Pd 

a 

P 

P 

> TJ 

P 

•H  P 

p 

00  VI 

44 

P 

*H  P 

Pd 

44 

p 

Vi  P 

/-V  Ti 

*0  P* 

TS  44 

c 

>^ 

p 

p 

P 

43  P 

•H 

P 0 

p 

P 

p 

w 0 

Vi  Pd 

3 

•H 

rH 

/-s 

V4 

44  dH 

Pd 

u d) 

p 

44 

44 

P 

P 

44 

p 

P 

P 

14-1  > 

dC  P 

P 

P 

P 

4J 

P 

0 

0 

3 

H 3 

0 

0 

P 

0 M 

p 

•H 

VI 

0 p 

0 

P 

Vi 

Pd 

0 P 

p 

4J 

0 M 

p 

SI 

S 

•H 

44 

V 

Vi 

0 

« P 

P 

a 

0 

fH 

Pd 

CM  44 

P 

P 

p 

A • 

VO  X 

0 

0 

dC 

iH 

rH 

V4 

fH 

0 

00  p 

1 

•H 

p 

p 

•H 

VI  ti 

cn 

P 

44 

•H 

M Vi 

44 

p 

44 

•H  P 

« 44 

rH 

P 

«d 

M P 

P 

X 

p 

0 P 

P 44 

rH 

44 

p m 

P 

•H 

Pu 

•H 

III.  ESTIMATION  OF  TANK  FAILURE  PRESSURE  j 

General  Discussion  and  Resulta  j 

I 

Estimates  were  made  of  the  minimum  Internal  pressure  loadings  i 

required  to  initiate  failure  of  cargo  tanks  of  three  representative 
vessel  designs:  an  offshore  barge,  an  Inland  barge,  and  a large 
tankshlp.  The  approach  used  to  develop  these  estimates  consisted  of 
elastic  structural  analyses  of  the  deck  and  sides  of  the  cargo  tanks, 
and,  where  appropriate,  the  bottom  and  bulkhead  structures.  The 
analysis  was  based  on  the  work  of  Clarkson  (1965)  on  plated  flat 
grillages,  augmented  by  the  orthotropic  plate  methods  of  Schade  (1940) 
and  Lekhnltskll  (1968). 

In  all  cases,  the  structural  resistance  of  the  plating  between 
the  longitudinal  and  transverse  framing  members  was  found  to  be 
considerably  higher  than  the  resistance  of  the  framing  members  them- 
selves, so  that  the  minimum  failure  pressure  loadings  were  defined 
by  the  stress  levels  developed  In  the  framing  members. 

The  Intent  of  this  analysis  was  to  determine,  using  simplified 
structural  models  of  the  cargo  tanks,  a range  of  pressure  levels  at 
which  typical  cargo  tank  configurations  would  be  expected  to  fall. 

These  pressure  levels  have  been  utilized  to  establish  the  overpressure 
hazards  for  the  range  of  cargo  venting  and  loading  rates  of  Interest 
to  this  overall  study  program.  It  was  not  the  intent  of  this  analysis 
to  conduct  a detailed  study  of  the  response  of  cargo  tank  structures 
to  the  various  static,  and  dynamic  loadings,  residual  stresses,  or 
thermal  effects  normally  considered  In  detailed  design  analysis.  Such 
additional  stresses  in  the  cargo  tank  structure  have  been  accoimted 
for  In  an  approximate  fashion  by  assuming  as  the  failure  criterion 
the  attainment  of  yield  stress  levels  (rather  than  ultimate  stress) 
in  the  tank  structure  due  to  the  internal  pressure  loading  alone.  The 
stresses  developed  due  to  the  other  loadings  were,  therefore,  assumed 
to  be  sufficient,  when  added  to  the  yield  stress,  to  produce  ultimate 
stress  levels  at  one  or  more  locations  In  the  highly  redundant  structure. 


The  accuracy  resulting  from  this  assimiptlon  and  from  the  elastic 
structural  analysis  Is  expected  to  be  about  + 30Z. 

For  the  purposes  of  this  study,  the  Coast  Guard  has  provided 
structural  data  on  the  cargo  tanks  of  the  three  representative 
vessels. 

On  the  basis  of  our  analyses  on  the  three  vessels.  It  has  been 
calculated  that  an  average  Internal  pressure  level  In  the  cargo  tanks 
of  about  six  pslg  will  be  sufficient  to  Initiate  failure  of  the  cargo 
tank  structure.  This  level  represents  an  average  of  about  el^t  pslg 
for  the  tankshlp,  six  pslg  for  the  offshore  barge,  and  about  four 
pslg  for  the  Inland  barge.  A summary  of  the  results  Is  given  In 
Table  III-l.  For  each  of  the  three  vessels  analyzed,  the  weakest  portion 
of  the  tank  structure  (the  probable  location  of  Initial  failure)  was 
determined  to  be  at  the  ends  of  the  transverse  frames  at  the  tank 
boundary  defined  by  the  deck-sldewall  intersection. 


B.  Structural  Analysis  Methods 


Many  ship  structures  are  In  the  form  of  plated  grillages.  In 
which  stiffening  beams  in  either  one  or  two  directions  (generally 
orthogonal)  are  attached  to  the  plating  to  provide  resistance  to 
lateral  loadings.  Such  structures  are  commonly  used  In  the  decks, 
sides,  bottoms,  and  bulkheads  of  ships;  and  also  find  applications 
In  civil  engineering  structures  and  In  aircraft  structures.  Grillages 
of  these  types  are  highly  redundant  structures,  so  that  they  cannot 
be  analyzed  by  static  equilibrium  methods  alone;  the  additional 
conditions  relating  to  the  compatibility  of  the  deflections  of  the 
structural  components  must  be  Imposed. 


For  most  grillage  structures  used  in  marine  cargo  tanks,  there 
are  generally  several  stiffening  beams  (longitudinal  and/or 
transverse  members)  so  that  the  overall  structure  has  a large  number 
of  redundancies.  Direct  solution  by  hand  of  the  stresses  and 
deflections,  using  indeterminate  analysis  methods.  Is  extremely 
cumbersome  for  grillages  of  any  significant  complexity,  with  the 


TABLE  III-l 


Cargo  Vessel 
Offshore  Barge 


Inland  Barge 


Tankshlp 


SUMMARY  OF  CARGO 

TANK  FAILURE 
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Structure 
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Pressure 
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Trans . Bulkhead 
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Trans . Bulkhead 

7 
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8 

Transverse  Beam 
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Based  on  average  of  upper  and  lower  limit  structural  models  of  tank 
structure  to  account  for  interior  trusswork.  See  text. 
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development  of  computer  methods,  however,  such  analyses  have  been 
carried  out  in  parametric  form  for  a range  of  structural  configurations. 
The  work  of  Clarkson  (1965)  has  included  study  of  flat  grillages 
under  uniform  lateral  pressure  (which  is  the  loading  of  principal 
concern  for  the  present  program)  for  several  cases  of  odd  numbers  of 
stiffening  beams,  up  to  orthogonal  sets  of  15  x 9 beams.  The  boundary 
conditions  studied  have  included  opposite  edges  either  fixed  or  simply 
supported.  The  plating,  which  in  most  ship  grillages  is  welded  to  the 
stiffening  members,  is  assumed  by  Clarkson  to  behave  as  an  effective 
flange  to  the  transverse  and  longitudinal  members,  so  that  the  structu- 
ral model  is  essentially  reduced  to  a grldwork  of  Intersecting  beams. 

One  of  the  difficulties  in  applying  Clarkson's  work  is  the 
relatively  limited  range  of  stiffness  ratios  studied.  Many  ship 
structures  of  current  design,  including  in  particular  bulkheads  in 
oil  tankers,  and  decks  and  sides  of  longitudinally-framed  tankers, 
consist  of  relatively  heavy  and  wide-spaced  transverse  members  and 
relatively  light  longitudinal  members.  In  addition,  some  bulkhead 
structures  are  typically  stiffened  in  one  direction  only. 

In  some  of  the  particular  grillage  structures  analyzed  in  this 
study,  it  was  found  that  the  Clarkson  data  could  be  extrapolated  to 
provide  reasonable  results;  in  other  cases,  the  extent  of  extrapolation 
needed  was  too  large  to  provide  meaningful  answers . An  alternative 
method  for  analysis  of  grillages  is  to  model  the  stiffened  structure 
as  an  orthotropic  plate,  in  which  the  stiffening  members  are  trans- 
formed into  a continuous  structure  having  an  equivalent  stiffness 
per  unit  width  in  each  of  the  two  orthogonal  directions.  Once  this 
is  done,  the  solutions  of  orthotropic  plate  analysis  can  be  applied, 
or  for  cases  in  which  the  stiffeners  and  their  spacing  are  the  same 
in  both  directions,  the  solution  of  the  flat  plate  analyses  can  be 
used. 

One  difficulty  with  the  orthotropic  plate  approach  is  that  only 
a limited  number  of  sets  of  boundary  conditions  have  been  formally 


solved.  The  solution  for  the  case  with  all  boundaries  fixed,  for 
example,  has  not  been  developed.  Estimates  for  the  fixed-boundary  case 
can  be  made,  however,  from  some  of  the  other  configurations  which  have 
been  studied,  both  theoretically  and  by  approximate  methods.  An 
orthotropic  plate  procedure  was  programmed  and  used  in  the  present 
study  to  develop  estimates  for  those  grillage  structures  for  which 
Clarkson's  data  sheets  could  not  be  used.  It  was  also  used  to  check 
the  results  for  those  cases  in  which  the  Clarkson  data  sheets  were 
applicable  either  directly  or  by  relatively  minor  extrapolations. 

C.  Offshore  Barge  Structural  Analysis 

1.  Dimensions 

The  offshore  barge  chosen  for  analysis  was  a Todd  605  design  of 
overall  dimension  376  feet  by  78  feet  by  31-3/4  feet.  The  largest 
cargo  tank  is  56  feet  long,  37  feet  wide,  and  31-3/4  feet  deep.  The 
deck,  sidewall,  bottom,  and  longitudinal  bulkhead  structures  of  the 
tanks  are  orthogonally-stiffened  with  relatively  small  and  close- 
spaced longitudinal  members,  and  with  relatively  heavy  transverse 
members  spaced  eight  feet  on  center.  These  transverse  members,  which 
are  notched  over  the  longitudinals  and  welded  directly  to  the  plating, 
are  part  of  a transverse  truss  structure  consisting  of  (1)  two 
stanchions  spaced  at  thirteen  foot  intervals  and  (2)  diagonal  members 
from  the  top  and  bottom  of  each  stanchion  to  the  center  height  of  the 
sidewall  and  the  longitudinal  bulkhead.  Such  transverse  truss 
structures  are  typical  of  many  barge  designs. 

The  transverse  bulkhead  structure  is  also  stiffened  in  both 
directions  with  closely-spaced  horizontal  members  and  heavy  widely- 
spaced  vertical  members. 

2.  Methods  of  Analysis  and  Results 

Each  of  the  major  structures  of  the  cargo  tank  was  analyzed  in 
two  ways:  with  the  grillage  methods  of  Clarkson  and  with  an  ortho- 
tropic plate  method.  The  presence  of  the  truss  systems  within  the 
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tank  complicates  the  modeling  of  the  deck,  bottom,  sidewall,  and 
longitudinal  bulkhead  structures  in  that  the  truss  contributes  signifi- 
cantly to  the  resistance  of  these  structures  to  lateral  pressure  loadings. 
An  accurate  determination  of  the  effect  of  the  truss  structure  in  this 
regard  would  require  a modeling  of  the  three-dimensional  configuration 
of  the  grillages  and  trusses  in  terms  of  one  of  several  available 
finite-element  computer  codes,  such  as  NASTRAN  or  STARDYNE.  Such  a 
procedure  is  clearly  well  beyond  the  scope  of  this  program.  An 
alternative  procedure  which  was  followed  here  was  to  attempt  to  bracket 
the  actual  structural  configuration  by  analyzing  simplified 
configurations.  These  simplified  models  Included  one  in  which  the 
out-of-plane  truss  members  were  ignored  (this  model  would  tend  to 
predict  much  lover  critical  pressures  than  the  actual  structure) ; 
another  in  which  the  transverse  section  of  the  plane  of  the  truss  was 
asstmied  to  be  Infinitely  stiff,  so  that  the  structure  analyzed,  by  an 
orthotropic  formulation,  consisted  of  a plate  reinforced  in  one 
direction  only,  and  having  a longitudinal  dimension  equal  to  the  truss 
spacing  (this  model  would  result  in  higher  critical  pressures  than  the 
actual  structure).  Still  another  model  analyzed  was  a grillage 
structure  whose  length  was  equal  to  the  full  tank  length,  but  whose 
width  extended  to  the  stanchion  only.  This  structural  approximation, 
which  was  based  on  the  assumption  that  the  truss  stanchions  defined 
a fixed  longitudinal  line,  would  also  result  in  prediction  of  higher 
critical  pressures  than  the  actual  structure  could  resist. 

As  an  example,  the  entire  deck  structure  of  the  barge,  when 
analyzed  without  the  contributions  of  the  interior  truss  members,  is 
calculated  to  develop  yield  stress  levels  in  both  the  longitudinal 
and  transverse  members  at  a uniform  lateral  pressure  of  about  two  psi. 

When  analyzed  with  the  transverse  sections  fixed,  yield  stress  levels 
would  develop  at  a uniform  lateral  pressure  of  about  nine  psi. 

Analysis  of  the  reduced  deck  structure,  * ssuming  a fixed  stanchion 
line,  results  in  a lateral  pressure  level  of  about  ten  psi  in  the 
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longitudinal  numbers  and  about  18  psl  In  the  transverse  members. 

A realistic  value  Is  obviously  above  two  and  below  nine,  and  has 
been  assumed  to  be  six  psl. 

Similar  approaches  used  for  the  bottom,  sidewall,  and  longitudinal 
bulkhead  structures  have  resulted  In  probable  values  for  the  critical 
lateral  pressure  equal  to  10  psl  or  more.  Calculations  for  the  trans- 
verse bulkhead  led  to  an  estimate  of  14  psl  for  the  critical 
pressure. 

For  completeness,  the  lateral  pressure  levels  necessary  to  develop 
yield  stress  In  the  Individual  areas  of  plating,  l.e.,  the  rectangular 
plate  segments  defined  by  the  grillage  members,  were  also  calculated. 
Assuming  fixed  boundary  conditions,  the  critical  pressures  were  In 
excess  of  12  psl. 

Thus,  the  critical  Internal  pressure  for  the  cargo  tanks  of  this 
offshore  barge  design  was  taken  to  be  six  psl,  corresponding  to  the 
development  of  yield  stress  In  the  barge  tankage  structures,  reached 
Initially  In  the  transverse  members  at  the  deck  sidewall  Intersection. 

D.  Inland  Barge  Structural  Analysis 
1.  Dimensions 

The  Inland  barge  design  studied  was  a slngle-skln  oil  barge, 

NBC  Hull  1654.  This  vessel  has  overall  dimensions  of  132  feet  by  50 
feet  by  10-1/4  feet,  with  a typical  cargo  tank  having  dimensions  of 
51  feet  by  25  feet  by  10-1/4  feet.  The  tank  structures  are  longi- 
tudinally stiffened,  with  transverse  frames,  also  In  the  form  of 
trusses,  spaced  at  approximately  7-1/2  foot  Intervals.  The  transverse 
members  are  welded  to  the  flanges  of  the  longitudinals , and  thus  are 
not  attached  directly  to  the  plating.  The  truss  system  Includes , In 
addition  to  the  top,  bottom,  and  sidewall  members,  three  stanchions 
spaced  at  6-1/4  feet  Intervals.  Diagonal  members  extend  from  the 
bottom  of  the  sides  and  from  the  bottom  of  the  middle  stanchion  to  the 
top  of  the  adjacent  stanchions,  thus  forming  an  Inverted  W across  the 


tank  section.  Both  the  longitudinal  and  transverse  bulkhead  structures 
are  stiffened  In  the  vertical  direction  only. 

2.  Methods  of  Analysis  and  Results 

Both  the  grillage  methods  of  Clarkson  and  the  orthotropic  plate 
approximations  were  used  to  analyze  the  various  structural  elements 
of  the  Inland  barge.  As  in  the  case  of  the  offshore  barge,  the 
analyses  of  the  deck,  bottom  and  sidewall  structures  were  based  on 
simplified  models;  the  first  model  neglecting  the  resistance  offered 
by  the  out-of-plane  truss  members,  another  model  based  on  the  assump- 
tion that  the  transverse  frame  was  infinitely  stxff,  and  a third  model, 
where  appropriate,  using  the  reduced  width  as  defined  by  the  stanchion 
spacing.  These  simplified  structural  models  can  reasonably  be 
considered  to  bracket  the  resistance  of  the  actual  grillage  and  truss 
structures . 

For  the  deck  structure,  the  analysis  of  the  entire  deck  without 
contributions  from  the  interior  truss  members  led  to  values  of  maximum 
lateral  pressure  of  about  one-half  of  a psi  for  initiation  of  yielding 
in  the  transverse  members,  and  about  one  psi  for  the  longitudinal 
members.  Assuming  that  the  transverse  members  (top  elements  of 
interior  truss  structure)  were  fixed,  yield  stress  levels  were  reached 
in  the  longitudinal  members  at  about  13  psi.  Finally,  the  model  based 
upon  the  width  of  deck  equal  to  the  spacing  of  the  stanchions  led  to 
critical  lateral  pressure  levels  of  about  eight  psi  for  the  transverse 
members  and  11  psi  for  the  longitudinals. 

A reasonable  mean  value  for  the  critical  lateral  pressure  for  the 
deck  structure  would  be  three  psi. 

Similar  sets  of  calculations  for  the  sidewall  and  bottom  structures 
indicated  that  the  critical  pressure  range  for  the  sidewall  was  3-1/2 
psi  for  yield  stress  in  the  transverse  members,  and  for  the  bottom,  a 
mean  value  of  about  four  psi  appears  reasonable. 
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The  bulkhead  structures  of  this  barge  design,  each  of  which  Is 
stiffened  In  the  vertical  direction  only,  were  analyzed  as  orthotropic 
plates,  and  led  to  critical  lateral  pressures  of  greater  than  seven 
psl.  The  plating  Itself  was  determined  to  require  In  the  worst  case 
a minimum  of  11  psl  to  reach  yield  conditions. 

The  rather  surprisingly  low  values  obtained  for  the  deck,  side, 
and  bottom  structures  of  this  barge  suggested  that  the  scantlings 
be  evaluated  for  adequacy  and  adherance  to  the  Rules  of  the  American 
Bureau  of  Shipping  (1967,  1971).  The  evaluation  Indicated  that  the  deck 
and  bottom  transverses  were  adequately  sized,  but  that  the  side  trans- 
verse was  undersized.  This  condition  was  brought  to  the  attention  of 
technical  staff  personnel  at  MMT  New  Orleans,  from  whom  the  scantlings 
had  been  obtained,  and  that  office  concurred  In  this  evaluation. 

In  view  of  this,  it  would  seem  reasonable  to  assume  that  the 
typical  Inland  barge  may  be  somewhat  more  resistant  to  internal 
pressure  loadings  than  Indicated  from  the  analyses  carried  out  for 
this  particular  barge  design.  Accordingly,  It  Is  suggested  that 
a value  of  the  critical  Internal  pressure  level  for  use  In  this 
study  equal  to  about  four  psl  would  be  appropriate. 

E.  Large  Tankshlp  Structural  Analysis 

1.  Dimensions 

The  tankshlp  chosen  for  analysis  was  a large  vessel  of  overall 
dimensions  1075  feet  by  143-1/2  feet  by  91  feet  deep,  with  a draft 
of  70  feet.  The  Inboard  deep  tank  analyzed  in  detail  was  198-1/2 
feet  long,  32  feet,  11  Inches  wide,  and  of  full  depth.  The  tank  was 
free-standing,  without  Internal  cross  members  or  trusses,  and  was 
longitudinally  framed  with  heavy,  widely-spaced  web  frames. 

2.  Methods  of  Analysis  and  Results 

The  deck  and  sidewalls  (which  In  the  c^e  of  the  inboard  tank 
also  serve  as  the  longitudinal  bulkhead  structures)  were  analyzed 
using  extrapolated  values  from  the  Clarkson  grillage  charts,  and  were 
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then  verified  using  the  orthotropic  plate  formulation.  Extrapolations 
to  the  Clarkson  charts  were  necessary  because  of  the  relatively  large 
web  frame  structures  used  in  the  tanks,  which  effectively  resulted  in 
a very  high  ratio  of  the  transverse  to  longitudinal  stiffness  of  both 
the  deck  and  sidewall  grillages. 

The  results  for  the  tank  ship  indicated  that  yield  stress  levels 
would  be  developed  at  the  upper  end  of  the  transverse  web  frame  of 
the  sidewall  under  an  internal  lateral  pressure  of  about  eight  psi, 
in  the  transverse  member  of  the  deck  structure  at  about  12  psi,  and 
in  the  longitudinals  of  the  deck  and  sidewall  structures  of  26  and  13 
psi,  respectively.  Yield  stress  levels  in  the  plating  itself  requires, 
as  in  the  cases  for  the  barge  tankage,  higher  lateral  pressures  than 
those  associated  with  the  development  of  yield  stress  in  the  grillage 
members . Assuming  fixed  boundaries , lateral  pressure  levels  of  24 
and  70  psi  would  be  needed  to  reach  yield  stress  in  the  sidewall  and 
deck  plates,  respectively.  On  the  basis  of  these  results,  it  would 
appear  that  a reasonable  assumption  for  the  value  of  the  critical 
internal  pressure  for  large  cargo  tank  structures  would  be  eight  psi. 
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IV.  TANK  PRESSURE  RISE  DURING  LOADING 


Factors  which  Contribute  to  Pressure  Rise 


Having  an  estimate  of  the  tank  failure  pressure,  we  seek  to  deter- 
mine the  maximum  pressure  reached  during  cargo  transfer  operations, 
and  when  it  is  reached.  As  illustrated  in  Figure  IV-1,  this  will  make 
it  possible  to  recognize  hazardous  cargo  transfer  situations.  The  ob- 
jective of  the  analysis  presented  is  to  predict  the  pressure  history 
inside  a tank  when  it  is  being  filled  with  a liquid.  The  pressure  rise 
results  from  the  finite  pressure  drop  of  the  vents  provided  to  discharge 
the  vapor/air  displaced  by  the  liquid  which  is  filling  up  the  tank. 


Estimated 
Range  of 
Tank  Failure 


Hazardous  Transfer  Operation 


Tank 

Pressure 


Safe  T ransfer  Operation 


Elapsed  Time  ■ 


FIGURE  IV-1  SCHEMATIC  OF  PRESSURE  HISTORY 
In  general,  (as  will  be  shown)  the  venting  capacity  is  adequate  to 
discharge  the  gas  mixture,  without  creating  back  pressures  in  excess 
of  a few  pslg.  However,  if  the  relief  valve  is  stuck  for  some  reason 
or  if  the  vent  area  is  reduced  to  blockage,  or  if  the  loading  rate  is 
excessive,  or  if  the  cargo  evaporates  vigorously  during  loading,  there 
results  a compression  of  the  vapor/air  mixture  in  the  tank  resulting  in 
hazardous  tank  pressures.  Some  of  these  factors  are  shown  in  Figure  IV-2. 


FIGURE  IV-2  SCHEMATIC  DIAGRAM  SHOWING  FACTORS  IMPORTANT  TO  PRESSURE  RISE 


The  phenomena  of  evaporation  and  mixing  of  vapor  with  air  during 
loading  affects  the  pressure  rise.  The  description  of  evaporation  Is 
In  general  exceedingly  complicated  and  demands  that  simplifying  assump- 
tions be  made.  Consider  the  bottom  loading  (of  liquid)  configuration 
shown  In  Figure  IV-3.  Just  before  the  liquid  Is  let  in,  the  tank  vol- 
ume Is  full  of  air.  As  the  liquid  comes  In,  It  seeks  to  establish  the 
equilibrium  vapor  pressure  next  to  the  liquid  surface.  Vapor  Is  gener- 
ated. If  the  vapor  is  heavier  than  air,  it  is  very  likely  that  there 
will  be  gas-air  stratification.  There  may  be  an  Intermediate  region 
where  there  will  be  a mixture  of  air  and  vapor  with  the  concentration 
varying  continuously.  On  the  other  hand,  if  the  vapor  generated  Is 
lighter  than  air,  mixing  caused  by  buoyancy  force  may  occur.  Thus,  the 
key  to  the  analysis  Is  how  the  mixing  process  Is  handled.  Molecular 
mlxln^’ls  negligible  compared  to  turbulent  mixing  as  a mechanism  for 
clearing  off  the  liquid  surface  of  vapor.  The  loading  pipe  outlet 
presumably  creates  enough  local  turbulence  to  justify  this  assumption. 

If  molecular  mixing  were  controlling,  there  would  be  negligible 
pressure  rise  and  evaporation  rate,  because  the  tank  dimension  greatly 
exceeds  the  distance  over  which  molecular  diffusion  can  penetrate  during 
the  fill  time. 

Turbulent  mixing  is  represented  by  a one-dimenionsal  Pick's  Law 

expression  for  vertical  diffusion.  The  coefficient  for  eddy  diffusion  Is 

4 2 * 

set  at  10  the  molecular  value:  (c  = 2 ft  /sec).  As  shown  In  a standard 

text  on  mass  transfer  (i.e. , Bird,  Stewart,  and  Llghtfoot,  p.  594  ff,  (I960)), 

one  can  derive  the  following  expression  for  the  mass  of  vapor  added  as  a 

function  of  time: 


M . S (1) 

vap  p V If  ^ ^ 

where  S Is  the  surface  area  of  liquid  and  the  equilibrium  partial 

pressure  of  vapor  (taken  at  the  liquid  cargo  temperature) . The  rate  of 
evaporation  Is  obtained  by  differentiation  of  Equation  (1)  with  respect 
to  time: 


This  is  quite  conservative  (more  mixing  than  actually  expected)  because 
the  eddy  transport  properties  are  commonly  taken  at  100  to  1000  times 
the  molecular  values. 

**Molecular  mixing  occurs  in  a stagnant  gas  and  thus  differs  from  turbulent 
mixing. 
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FIGURE  IV-3  SOME  PHYSICAL  PROCESSES  WHICH 
AFFECT  THE  EVAPORATION  RATE 
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(2) 


The  analysis  was  carried  out  omitting  the  correction  term  2t(dp/dt)/p, 
which  ranges  between  .002  and  .02  for  the  parameter  values  of  interest. 
Furthermore,  the  pressure  p in  the  denominator  of  Equation  (2)  was  taken 
to  be  a constant  equal  to  the  tank  pressure  when  loading  is  complete, 
p(tfiii).  This  is  consistent  with  omitting  the  correction  term 
(dp/dt  = 0 implies  p = constant).  Essentially,  the  net  effect  of  these  two 
changes  is  to  omit  the  effect  on  evaporation  of  the  Increase  in  ullage  pres- 
sure during  the  cargo  loading  process. 

* 

Since  Equation  (1)  suggests  '''  P > taking  p = will  under- 
estimate the  mass  of  vapor  added  since  = p(t).  However, 

this  simplification  is  justifiable  in  view  of  the  other  coarse  assumptions 
Imposed  on  the  evaporation  analysis  (see  below) . The  actual  variation  of  p 
which  we  have  omitted  in  Equation  (1)  is  much  less  than  the  uncertainty 

in  the  diffusion  coefficient  e.  The  expression  actually  adopted  in  the 
analysis  is  based  on  evaporation  into  a tall  uncovered  tank: 


K 

vap 


p p S 
vap  vap 


(3) 


In  the  calculations  of  tank  pressure,  it  will  be  convenient  to  use  the 

ratio  K of  volumetric  evaporation  rate  (M  /p  ) to  the  volumetric  loading 
. vap  vap 

rate  Q.  The  ratio  K must  be  referenced  to  a particular  time  since  M 

-1/2 

decreases  as  t . For  convenience,  we  take  the  reference  time  as  the 


fill  time,  t 


fill' 


K = 


/ EVAP.  RATE  \ 
VlOADING  RATEy  ^ 


fill 
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A sample  value  of  K can  be  obtained  as  follows:  substituting  e = 2 ft  /sec, 

Q - 5 ft^/sec,  - A hr,  S * 500  ft  , p^^^  = 5 psi,  and  * 6 psi, 

we  obtain 


K = .55 


t 


1 


Values  of  K ranging  from  0.1  to  0.7  were  considered  in  the  analysis  in  order 
to  explore  larger  values  of  the  liquid  surface  area  S and  the  diffusion 
coefficient  c. 


Substituting  Equation  (4)  into  Equation  (3),  and  assuming  p = p 


we  obtain  the  working  expression  for  evaporation  rate: 

iap  ■>'  “a  ^ 


vap’ 


(5) 


Additional  complications  arise  (a)  if  the  vapor  is  hotter  or  colder 
than  initial  air,  (b)  if  there  is  a cooling  effect  of  evaporation  due 
to  latent  heat,  or  (c)  if  the  actual  exposed  liquid  surface  area  varies 
due  to  a jet  from  the  fill  pipe  or  spreading  across  the  tank  floor. 

Here  we  neglect  all  of  these  complications  by  assuming  that: 

• Cargo  loading  temperature  is  close  to  ambient; 

• Liquid  surface  stays  at  the  cargo  loading  temperature  (which 
is  known) ; and 

• Constant  exposed  liquid  surface  area  equal  to  tank  floor  area. 

C.  Mathematical  Model  of  the  Venting  Process 

Mathematical  models  have  been  developed  which  describe  the  pressure 
as  a function  of  time  for  given  cargo  transfer  conditions.  The  model 
presented  below  Includes  the  effects  of  gas  compressibility,  vent  system 
friction  and  restrictions,  cargo  evaporation,  loading  rate,  and  heat  trans 
fer  from  the  ambient  to  the  gas  in  the  vent  system. 

1.  Assumptions 

(a)  Gas  properties:  The  entire  analysis  is  performed  assuming 
that  the  vapor/air  mixture  being  vented  has  the  molecular  weight, 
specific  heat,  and  viscosity  of  pure  air.  In  the  case  of  high  vapor 
pressure  liquids,  there  will  be  substantial  evaporation  into  the  tank, 
and  if  the  vapor  differs  substantially  from  air,  this  assumption  will 
not  be  strictly  valid. 

(b)  Liquid  loading  rate  is  constant. 

(c)  Tank  volume  is  fixed  - l.e.,  we  neglect  tank  expansion. 

(d)  The  friction  coefficient  for  the  venting  system  is  a constant 
(that  is,  the  flow  through  the  vent  is  in  the  high  Reynolds  Number  - 
turbulent  flow  regime). 


t 


(e)  Vent  flow  is  isothermal:  Heat  exchange  occurs  between  the 
ambient  and  the  gas  flowing  in  the  vent  system  to  the  extent  that  the 


gas  temperature  variations  are  negligible.  An  alternate  model  based 
on  adiabatic  flow  is  presented  in  Appendix  A. 


(f)  Evaporation  rate  is  limited  by  turbulent  vertical  mixing 
in  the  tank,  so  that  the  mass  addition  rate  is  proportional  to  the 
inverse  square  root  of  elapsed  time , 

(g)  Initial  conditions  within  the  tank  are  the  same  as  the  ambient 
conditions. 

(h)  The  mass  of  vapor  is  assumed  small  compared  to  the  liquid 
cargo  mass,  even  at  earliest  stages  of  loading. 

2.  Analysis 

First  a mass  balance  equation  is  written  for  the  vapor  mass  in  the  tank 
as  a function  of  the  outflow  rate  and  the  evaporation  rate.  Then  from  the 
vent  characteristics,  the  outflow  rate  is  expressed  as  a function  of  the 
tank  pressure.  The  evaporation  rate.  Equation  (5),  is  also  used  in  the 
mass  balance.  The  resulting  equation  for  tank  pressure  is  complicated 
and  must  be  solved  simultaneously. 

Taking  a mass  balance  on  the  vapor/alr  mixture,  we  obtain: 
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Equation  (6)  is  the  fundamental  equation  determining  the  pressure 
rise.  Each  of  the  three  terms  will  be  expressed  in  terms  of  tank  pres- 

I 

sure.  The  left-hand  side  is  developed  as  follows.  For  a compressible  j 

gas,  the  mass  of  the  air-gas  mixture  is 


For  constant  loading  rate,  V/V^  = 1 - isothermal 

gas,  p/p  » p/p,;  therefore.  Equation  (7)  becomes 

a a 

M-  p^VT(P/Pa)(l  - t/t^^^^) 


*3Q 


(7) 


(8) 


Next,  the  venting  rate  is  to  be  related  to  the  pressure  rise. 

The  relationship  between  pressures  on  the  upstream  and  downstream  side 
of  the  vent  and  the  mass  flow  rate  can  be  derived  for  isothermal  flow  from 
Shapiro, eq.  (6.42)  on  page  182: 


2 2 

P - Pa 


■in 


RT  [2  Jin  (P/Pg)  + -g] 
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where  R is  the  gas  constant.  Values  of  the  Fanning  friction  factor 
for  fully  developed  turbulent  flow  under  typical  venting  conditions 
range  from  0.004  to  0.006  with  a typical  value  of  about  0.005. 

Substituting  P/p  = p/Pa  and  solving  for  M^,  we  obtain 


^ ^ (P^  - Pa^)  _ A 

^ (2RT)^^^[Jln(p/p3)  + 2f  L/D]^/^ 

In  order  to  complete  the  formulation,  we  recall  from  Section  IV-A  that 


M =KM(t  /t)- 
vap  c fill 


where  K is  defined  on  page  28. 

Substituting  Equations  (5),  (8(,  and  (10)  into  Equation  (6), 
we  obtain  an  expression  for  tank  pressure  ratio  p/p^: 


^a'^T  dt  f^P^'Pa^^^  " ‘^''‘^fill^^ 


[(p/pa^^  - 1]^^^  A Pa 
(RT)^^^[2Jln(p/pa)  + 4f  L/D]^^^ 


Differentiating  the  left  hand  side  and  dividing  by  m^  , we  can  derive 


t (1  - t/t  ) = p/p  + K 

*^flll  fill  dt  ^ 


(M' 


ii 


2 _ 
f(p/pa)  - 11  /RT 

(Q/A)[2!ln(p/Pa)  + L/D]‘ 


* R is  defined  as  follows:  R • .00232/ (molecular  weight)  in  units 

p8i2-ft^-sec2-lbm-2-®R-l.  The  acceleration  due  to  gravity  is  included 
in  the  gas  constant  to  convert  Ibm/in^  to  psl. 


Equation  (12)  is  subject  to  the  initial  condition  p/p^  = 1 at  t - 0 
and  can  be  numerically  integrated.  The  initial  problem  of  K/ /t  going  to  “> 
as  t "*■  0 was  overcome  by  a numerical  technique.  The  key  parameters  in  Equa- 
tion (12)  are  Q/A,  K,  and  4f  L/D.  Values  of  these  parameters  were 

selected  from  Table  T.I-3  (page  13). 

D.  Results  for  Non-Volatile  Cargoes  (K  = 0) 

In  order  to  compare  the  pressure  time  history  expected  for  various 
vent  systems,  cargoes,  and  loading  rates,  ranges  of  values  of  the  different 
cargo  transfer  parameters  were  chosen  from  Table  II-3  as  follows: 

4f  L/D  ratio  values  = 20,  60,  and  200; 

Q/A  values  = 10,  30,  and  100  ft /sec. 

It  was  found  from  the  results  that,  for  "typical"  Q/A  and  4f  ^ 
values,  the  predicted  pressure  rise  is  very  small.  The  predicted 
pressure  rise  obtained  under  the  conditions  4f  L/D  = 10  and 
Q/A  = 100  ft/sec  is  shown  in  Figures  IV-4  and  IV-5.  In  these 
figures,  the  predicted  tank  pressure  rise  is  plotted  against  the  fraction  of 
total  fill  time.  Also  indicated  in  the  figures  is  the  pressure  range  in 
which  the  tank  is  expected  to  fail.  In  Figure  IV-4,  the  pressure  is 
plotted  on  a log  scale  to  bring  out  the  differences  in  the  predicted 
pressure  values  clearly.  In  Figure  IV-5,  linear  scale  plots  are  used 
so  that  the  timewise  development  of  the  pressure  is  shown  more  clearly. 

In  addition,  the  relative  magnitudes  of  the  highest  tank  pressure 
achieved  and  the  expected  tank  failure  pressure  range  are  shown. 

Table  IV-1  shows  the  maximum  pressure  as  obtained  from  incompres- 
sible flow  model*and  the  compressible  adiabatic  flow  model  for  different 
characteristic  Q/A  velocities  and  4f  L/D  values.  It  is  seen  that  for 
several  of  the  cases  considered  with  compressible  flow,  the  Mach  Number 
at  the  exit  section  becomes  1,  resulting  in  choked  flow  conditions. 

When  such  a situation  results,  the  calculation  procedure  programmed  on 
the  computer  leads  to  physically  Impossible  solutions  (because  in  the 
adiabatic  calculations,  provision  is  not  made  for  the  occurrence  of 
choked  flow).  These  situations  are  Indicated  in  Table  IV-1  with  "X" 
marks . 

’*The  Incompressible  flow  model  is  the  same  as  equation  (12)  for 
P/Pa=l. 
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It  is  seen  from  Figure  IV-4  that  there  is  not  a lot  of  difference 
in  the  maximum  pressure  values  from  incompressible  vs.  adiabatic  com- 
pressible flow.  This  is  understandable  because  even  under  the  condition 
of  Q/A  = 100  ft/sec,  the  flow  Mach  Number  is  only  about  0.1.  This  indicates 
that  compressibility  effects  are  small. 

The  second  important  point  noticed  from  the  calculations  is  that 
the  maximum  pressures  are  far  less  than  the  failure  pressure  of  the 
tank.  This  indicates  that  the  venting  system  (even  at  the  high  gas 
flow  velocities  of  100  ft/sec)  is  capable  of  venting  all  of  the  gas 
that  is  being  expelled,  with  pressure  build-up  not  more  than  a few 
psi.  It  is  emphasized  that  the  results  hold  good  only  for  gas  venting 
and  do  not  extend  to  liquid  venting. 

Results  shown  in  Table  IV-1  Indicate  that  when  either  flow 

velocity  is  high  or  friction  is  high,  the  back  pressure  (tank  pressure) 

becomes  high  and  the  differences  using  incompressible  and  compressible 

• ^ 

become  considerable.  For  example,  for  Q/A  = 100  ft/sec  and  L/D  = 3000, 
the  peak  pressure  predicted  by  incompressible  methodology  gives 
= 4.85  psig,  whereas  with  the  compressible  methodology,  = 10.5  psig 

(which  is  2.17  times  that  given  by  incompressible  theory  on  gage  pressure 
basis) . This  indicates  that  compressibility  of  gas  becomes  a very 
important  issue  when  the  frictional  resistance  in  the  venting  valve-pipe 
system  becomes  large.  It  is,  however,  comforting  to  note  that  for 
typical  shiptank  venting  systems,  (L/D)  is  of  the  order  of  150  to  500. 

In  such  cases,  the  peak  pressure  can  be  predicted  very  easily  using  the 
incompressible  equation. 

From  the  analyses  and  results,  we  conclude  that  for  non-volatile 
cargoes  (K  = 0) , 

• The  peak  tank  pressures  resulting  during  the  venting  of  gas  are 
small  compared  to  the  tank  failure  pressures. 

• The  peak  pressures  can  be  calculated  by  a simple  incompressible 
analysis  when  the  gas  flow  rates  in  the  vent  system  are  small 
compared  to  the  velocity  of  sound  in  the  ambient  atmosphere. 


Incompressible,  No  Evaporation 


max 


(PSIG) 


Compressible,  No  Evaporation 


• High  friction  (caused  by  a constricted  vent  or  extremely  long 
length  of  venting  pipes)  tends  to  make  the  peak  pressures  high. 
Under  these  conditions,  the  value  obtained  from  simple  incom- 
pressible analysis  is  inaccurate. 

E.  Results  for  Venting  with  Cargo  Evaporation 

1.  Magnitude  of  the  Pressure  Rise 

For  normal  vent  systems,  cargo  volatility,  and  loading  rates,  it 
appears  that  there  is  about  a safety  factor  of  three  in  tank  pressure  rise. 
For  example,  we  have  solved  Equation  (12)  for  the  case  of 

• Vent  restriction  4f  L/D  = 20; 

Loa.d.ln|  Rate  . ' , 

Vent  Area 


Evap  Rate 
Loading  Rate 


= K = 0.5 


For  this  case  the  pressure  Increases  rapidly  to  about  2 psig  during  the 
first  4%  of  the  loading  time  (approximately  10  minutes).  Then  as  the 
evaporation  rate  falls  off,  the  tank  pressure  drops  to  less  than  1 psig. 

In  practice  the  relief  valve  would  close  as  the  pressure  dropped  and  re- 
open only  for  pressures  greater  than  the  set  point  (2  psig  is  common). 
Comparison  of  this  result  to  the  results  in  Section  D for  no  evaporation 
indicates  that  evaporation  has  the  effect  of  increasing  the  maximum 
pressure  early  in  the  loading  period  by  about  1 psig.  This  does  not  pre- 
sent a hazard,  but  again  only  holds  for  normal  venting  conditions.  Let 
us  examine  the  individual  factors  which  can  lead  to  higher  pressure  rise 
during  abnormal  cargo  transfer  situations. 

2.  Effect  of  Vent  Restriction 

Normally,  the  vent  system  restriction  (with  valves  open)  corresponds 
to  an  effective  L/D  of  500  for  manifold  systems.  For  f » .005  this  cor- 
responds to  4f/L/D  = 10.  Again,  this  normal  extent  of  vent  restriction 
gives  rise  to  tank  pressures  of  less  than  2 psig.  Now  if  L/D  is  Increased 
a factor  of  6,  the  peak  pressure  goes  up  to  3.3  psig  and  is  reached  at 
about  14%  of  the  total  fill  time,  as  shown  in  Figure  IV-6.  If  the  ef- 
fective L/D  is  Increased  a factor  of  20,  say  by  a blockage  or  sticky  PV 


Fraction  of  Fill  Time 


FIGURE  IV-6 


EFFECT  OF  VENT  RESTRICTION 


I 
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valve,  then  the  maximum  pressure  Is  over  7 pslg  and  Is  stabilized  at 
about  40%  of  the  total  fill  time.  Clearly,  tank  failure  Is  risked  only 

for  restrictions  having  an  effective  L/D  on  the  order  of  6,000  which  Is 

more  than  ten  times  normal  for  a manifold  vent  system. 

3.  Effect  of  Loading  Rate  and  Vent  Capacity 

The  key  parameter  here  Is  the  ratio  of  loading  rate  to  vent  cross 

sectional  area,  Q/A,  which  has  the  units  of  velocity  and  Is  normally 

10-30  ft/sec.  If  for  some  reason  the  loading  rate  Is  Increased  by  a 
factor  of  3.3  (from  30  ft/sec  to  100  ft/sec),  then  the  peak  pressure 
would  go  up  from  2 to  6.5  pslg,  as  shown  In  Figure  IV-7.  This  Is  on 
the  border  of  the  hazardous  range.  If  the  loading  rate  Is  ten  times 
normal,  or  If  the  vent  pipe  size  Is  one-tenth  the  usual  area,  then 
the  tank  pressure  goes  up  steadily  and  can  cause  tank  failure  after  less 
than  10%  of  the  fill  time  (24  minutes).  The  Q/A  and  L/D  combinations 
leading  to  hazardous  tank  pressures  are  to  be  found  from  Table  IV-2. 

4.  Effect  of  Cargo  Volatility 

The  effect  of  cargo  volatility  manifests  Itself  In  the  first  one- 
fifth  of  the  cargo  transfer  operation,  as  shown  In  Figure  IV-8.  Using 
Identical  loading  rate  and  vent  system,  a cargo  with  a high  vapor 
pressure  (K  = 0.7)  gives  rise  to  a peak  tank  pressure  of  2.6  pslg 
compared  to  the  value  of  1.0  pslg  reached  by  a less  volatile  cargo 
(K  = 0.3).  Note  that  these  tank  pressures  are  much  less  than  the 
corresponding  equilibrium  vapor  pressures.  Implying  that  diffusion 
markedly  limits  the  extent  of  vapor  addition.*  The  tank  pressures  are 
well  below  the  hazardous  regime  above  6 pslg,  so  that  It  appears  that 
even  an  extremely  volatile  cargo  (p  = 14,7  psl  at  70°F)  would  not 
cause  a problem  unless  Q/A  or  L/D  Is  excessive. 


*Thls  comment  Is  subject  to  the  assumptions  of  the  evaporation  model. 
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FIGURE  IV-7  EFFECT  OF  LOADING  RATE 
AND  VENT  CAPACITY 


L/D 


L5 


20 


25 


30 


Q/A  (FT/SEC) 


40 


600 

Tank  pressure 

less 

than 

/ 

300 

3 psl 

/ 

3.0 

IJOO 

/ 

3.3 

1200 

/ 

2.8 

3.6 

1600 

/ 

3.1 

4.1 

2000 

/ 

3.0 

3.5 

4.6 

2400 

/ 

3.2 

3.8 

5.0 

32JO 

/ 3.0 

3.6 

4.3 

5.9 

4000 

_/ 

3.3 

4.0 

4.8 

6000 

2.2 

3.0 

4.0 

4.9 

6.0 

OJOO 

2.5 

3.5 

4.6 

5.8 

10000 

2.7 

3.8 

5.1 

/ 

12000 

2.9 

4.2 

5.6 

/ 

1'  JOO 

3.3 

4.8 

20000 

3.6 

5.4 

/ 

24000 

3.9 

6.0 

/ 

/ 

320)0 

4.5 

/ 

4J000 

5.1 

/ 

50 


5.1 


5.8 


60 


4.1  4.9 


4.4  5.4 


80  100  150  20 


2.7  3.2  4.2  5.3 


3.2  3.8  5.1 


4.4  5.9 


Tank  pressure  greater 
than  6 psl  - hazardous 
for  most  vessels 


TABLE  IV- 2 

Maximum  Tank  Pressures  (PSIG)  for 
Given  Values  of  Vent  System 
L/D  and  Q/A  (Based  on  f = .005) 
K • 0.7 
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F.  Procedure  for  Evaluating  Overpressure  Hazard  in  an  Individual  Tank 

2 

1.  Determine  the  vent  area  A (ft  ) from  the  cross-sectional  area  of  the 
vent  pipe(s)  leading  from  the  tank. 

. 2 

2.  Determine  the  loading  rate  Q (ft  /sec)  for  the  tank  in  question  from 
the  total  pumping  rate  and  the  number  of  tanks  to  be  loaded  simul- 
taneously. 

3.  Calculate  the  quotient  Q/A  (ft/sec)  which  is  the  characteristic 
venting  velocity. 

4.  Determine  the  effective  length-to-diameter  ratio  L/D  of  the  vent 
system  using  Table  II-2  and  the  procedure  illustrated  in  Figure  II-6. 

5.  Use  Table  IV-2  to  find  the  maximum  tank  pressure  to  be  expected 
during  cargo  transfer,  using  the  values  of  L/D  and  Q/A  obtained 
in  steps  3 and  4. 

■ 6.  Compare  the  expected  tank  pressure  to  the  failure  pressure  which 

may  be  obtained  (1)  from  a detailed  structural  analysis  of  the 
vessel  in  question  (see  Chapter  III)  or  (2)  assumed  to  be  6 psig 

, in  lieu  of  structural  data. 
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V.  TANK  PRESSURE  HISTORY  DURING  OVERFILL 
^ A.  Problem  Formulation 

The  objective  is  to  predict  the  tank  pressure  history  when  the  tank 
is  overfilled  while  being  loaded  with  liquid  cargo.  The  physical  system 
to  be  modeled  is  shown  in  Figure  V-1.  At  the  point  at  which  the  tank  is 
full,  the  pumping  process  will  begin  to  displace  liquid  into  the  vent 
system.  This  causes  a pressure  drop  due  to  friction  and  due  to  the 
vertical  head  of  the  liquid  in  the  vent  piping.  The  system  attempts  to 
relieve  this  pressure  drop  in  two  ways:  (a)  by  compressing  the  liquid, 
and  (b)  by  expansion  of  tank  walls. 

In  the  present  analysis,  the  tank  pressure  is  defined  as  the  pressure 
in  the  liquid  at  the  top  part  of  the  tank.  As  shown  in  Figure  V-2,  the 
pressure  at  the  top  of  the  tank  is  less  than  the  pressure  at  the  bottom 
of  the  tank  because  of  the  liquid  head. 

Additional  Pressure  Drop 
From  Fittings 


t = 0 

Cargo  Fills  Tank 


t = 1 - 1 0 sec. 
Cargo  Enters  Vent 


t = 20  - 200  sec. 
Cargo  Fills 
Vent  System 


FIGURE  V-1  SCHEMATIC  OF  THE  PHYSICAL  SYSTEM  FOR 
MODELING  THE  TANK  PRESSURE  RISE 
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Pressure 


FIGURE  V-2  VERTICAL  PRESSURE  GRADIENT 
DUE  TO  HYDROSTATIC  HEAD 


B.  Analysis 

We  first  write  a mass  balance  for  the  tank-liquid  system.  Then  the 
rate  of  mass  outflow  (exit  velocity)  and  the  rate  of  mass  accumulation 
within  the  tank  are  related  to  the  tank  pressure.  This  results  in  an 
equation  describing  the  time  rate  of  increase  of  pressure,  which  is 
solved  numerically  to  obtain  tank  pressure  as  a function  of  time.  These 
steps  are  shown  below. 


dt 


pAU 


Rate  of  accumulation 
of  mass  within 
the  tank 


Inlet  mass 
flow  rate 
of  liquid 


Outflow  rate 
of  liquid 
from  the  tank 


Here  we  take  U to  be  the  average  liquid  velocity  in  the  vent  system  at 
any  instant.  The  rate  of  accumulation  of  mass  is  related  to  the  rate 
of  pressure  rise  through  equation  (18),  which  is  derived  as  follows: 

M = pV^ 

where  is  the  tank  volume. 

Differentiating  equation  (14)  with  respect  to  time,  we  get 

dM,  p^!l+ 

dt  dt  Tdt 

For  constant  temperature,  we  have 


p = p(p) 


Therefore, 


dt  dt  L dp  ^ dp  J ^dt  LVm  dp  P dpj 


Using  the  definitions  of  compressibility  and  the  volumetric  expansion 
coefficient  (see  nomenclature),  we  have: 

f 

Substituting  equation  (18)  into  equation  (13) , we  get 

d£-  1 

dt (6  + k)  I M J 

I 

The  mass  in  the  tank  M is  given  by 


t 

M = M(0)  + / (M^  - pAU)  dt 


Hence,  equation  (19) can  be  written,  after  substituting  for  M from  equation  (20) 
and  dividing  both  numerator  and  demonlnator  by  pA,  as 
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dp 1 

dt  “ (B  + k) 


(Q/A)  - U 


(V  /A)+  (Q/A)  t - / U dt 


(21) 


Before  we  can  solve  equation  (21 ),  we  will  need  to  relate  the  outflow  velocity 
U with  the  tank  pressure  p. 

The  above  equation  is  the  fundamental  equation  for  estimating  the  pressure 
in  the  tank  as  a function  of  time.  Depending  on  the  assumptions  for 
U (p),  B»  and  k,  the  rate  of  rise  of  pressure  will  be  different. 

The  following  assumptions  and  effects  are  reasonable  for  typical 
cargoes  and  vent  system  designs: 

1.  Constant  mass  inflow  rate  ft-^.  The  change  in  liquid 
density  due  to  compression  has  a negligible  effect  on 
the  loading  rate  (for  example,  the  change  in  the  density 
of  Hexane  is  0.04%  for  every  atmosphere  change  in  pressure). 


2.  Constant  tc  value. 


3.  Expansion  of  the  tank  due  to  the  internal  pressure  taken  into 
account  by  a linear  expression  (constant  value  of  B) . 

4.  The  flow  in  the  vent  pipe  is  fully  turbulent  (Re  > 3000) 
and  therefore  the  friction  factor  is  a constant. 

5.  Increased  pressure  drop  occurs  as  the  length  of  liquid 
wetting  the  pipe  Increases. 

6.  A sudden  increase  in  pressure  drop  occurs  due  to  the  flow 
resistance  through  the  valve  and/or  flame  arrestor. 


The  pressure  drop  across  the  vent  system  will  be  due  to  both  ver- 
tical head  of  liquid  and  frictional  drop.  Therefore,  the  vent  system 


must  be  described  before  deriving  Ap.  The  vent  system  is  assumed  to 
consist  of  a vertical  pipe  section  of  height  H followed  by  a horizontal 
(or  near  horizontal)  section  ending  in  a P-V  valve  and  a flame  arrestor, 
as  shown  in  Figure  V-3.  The  total  geometric  length  of  the  pipe  equals 


L 


1" 


The  total  effective  length,  including  valves  and  fittings,  is  L. 


^ I 
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FIGURE  V-3  SCHEMATIC  DIAGRAM  OF  THE  VENT  SYSTEM 


The  tank  pressure  is  related  to  mass  outflow  velocity  through  the  vent  system 
by  the  expression 

p = p^  + p gx+-|  (4f)  Q P (22) 

where  x is  the  vertical  column  of  the  liquid  in  the  vent  pipe  (x  ^ H) 
and  y = total  length  of  liquid  wetted  pipe  at  any  instant  (y  < L]^)  • 

In  using  Equation  (22),  the  reader  is  cautioned  to  divide  the  second  and 

2 

third  term  by  144  in  order  to  convert  Ibm/ft  to  psi  which  appears  through- 
out the  text. 

When  y = liquid  has  completely  filled  up  the  pipe  and  is  venting 
out  to  the  atmosphere.  Under  this  situation,  the  length  y which  appears 
in  the  last  term  of  equation  (22)  becomes  L,  the  effective  length  of 
the  vent  passage,  including  the  equivalent  length  of  valves  and  fittings. 

Thus,  the  definitions  needed  to  supplement  equation  (22)  are  as  follows: 

jy  for  V < H 

(23) 

H for  y > H 


y = 


dt  for  y < L 


1 


L otherwise 


(24) 


Equations  (21) and  (22)  form  a coupled  set  of  two  nonlinear  equations 
for  two  unknowns,  U and  p.  It  is  also  necessary  to  use  equations  (23^ 
and  (24)  in  the  solution. 

The  initial  conditions  are  given  by 

p = p^  at  t = 0 where  p^  is  given  by  solving  equation  (12)  at 
t = ^flll  for  gas  venting. 

U = Uo  at  t = 0 

The  iterative  solution  procedure  for  obtaining  p and  U as  functions  of 
time  is  shown  in  Appendix  C.  The  results  for  specific  tank  conditions 
are  discussed  in  a later  section.  The  initial  condition  for  vent  exit 
velocity  is  derived  in  Appendix  D.  The  key  parameters  in  the  solution 
are  k,  H,  Q/A,  and  4f  L/D. 

C . Determination  of  Parameters 

Suppose  a tank  17  ft  wide  by  40  ft  long  by  90  ft  deep  is  being 
filled  with  hexane.  The  vent  pipe  is  4 inches  in  diameter  and  50  ft 
long  and  has  a vertical  section  which  is  10  ft  in  length.  The  liquid 
filling  rate  is  such  that  the  entire  tank  fills  up  in  four  hours.  The 
equivalent  (L/D)  ratio  for  friction  through  the  valve  is  100.  The 
data  for  a chemical  tanker  needed  to  derive  the  coefficients  for  equations 
(21)  - (25)  are  as  follows: 

3 

• Volume  of  tank  = V.j,=  61200  ft 

• Tank  filling  time  = V.j,/Q  = 4 hours 

• 3 

• Inflow  rate  (volumetric)  = Q = 4.25  ft  /s 

• Total  length  of  vent  pipe  = 1^=  50  ft 

• Diameter  of  vent  pipe  = D = 1/3  ft 

• Equivalent  (L/D)  for  valve  = (L/D)^a]^^g  = 100 

-3 

• Volumetric  expansion  coefficient  for  the  tank  » 3 • 1.7  x 10  /atm 


(see  Appendix  E) 

• Coefficient  of  friction  for  flow  in  vent  pipe  *4f  - 0.01 


• Liquid  used  in  the  tank  = Hexane 

• Density  at  ambient  pressure  and  temperature  = =41.12  Ibm/ft 

• Compressibility  = k = 2.75  x 10  psi 

(see  Handbook  of  Physics  and  Chemistry,  page  209,  Table  165) 

The  initial  flow  velocity  (at  zero  time)  of  liquid  into  the  vent 
pipe  is  calculated  using  the  procedure  illustrated  in  Appendix  D . The 
value  obtained  for  the  above  values  of  the  parameters  is 

U =1.6  ft/sec. 


D.  Results  and  Discussions 


1.  Hazardous  Loading  Conditions 

The  liquid  efflux  velocity  U which  corresponds  to  a hazardous  6 
psig  tank  pressure  can  be  obtained  from  Equation  (22) : 

(p  - p - p g 

1 ^ 

^ P 4f  L/D 

Substituting  4f  L/D  = 10  and  H = 8 ft  for  a normal  manifolded  vent  system, 
we  obtain  U = 9 ft/sec.  If  Q/A  is  less  than  approximately  9 ft/sec,  pres- 
sures approaching  those  necessary  to  rupture  the  tank  would  not  be  expected. 
Liquid  cargo  can  pass  through  the  vent  system  indefinitely  only  if  the  vent 
capacity  and  loading  rate  are  such  that  U is  well  below  this  value. 

For  higher  loading  rates  or  smaller  vent  capacities,  tank  failure 
appears  to  be  inevitable.  The  pertinent  question  is  how  much  time  does 
the  operator  have  to  react  to  an  overfill  condition.  We  have  conducted  a 
numerical  analysis  for  two  cases  in  this  failure  regime: 

o Loading  rate/vent  area  = Q/A  = 49  ft /sec;  and 

o Loading  rate/vent  area  = Q/A  = 98  ft/sec. 

The  results  of  the  numerical  integration  of  Equations  (21)  and  (22)  are 
shown  in  Figures  V-4  and  V-5.  The  tank  (top)  pressure  is  plotted  as  a 
function  of  dimensionless  time.  Figure  V-4  is  for  a four-hour  fill  time 
and  V-5  for  a two-hour  fill  time.  Results  for  three  values  of  the  coef- 
ficient 4f  L/D  appear  on  each  Figure:  10,  100,  and  “ (blocked  vent). 

From  Table  II-3  it  can  be  seen  that  typical  vent  systems  are  in  the  range 
3 < 4f  L/D  < 10. 

2.  Effect  of  Vent  L/D,  Cargo  Viscosity,  and  Vertical  Runs 

It  is  seen  from  Figures  V-4  and  V-5  that  most  of  the  pressure  rise 
is  due  to  the  friction  generated  by  flow  in  the  vent.  The  effect  of  hydro- 
static head  is  negligible  (hydrostatic  pressure  due  to  10  feet  of  vertical 

pipe  with  Hexane  is  approximately  3 psig). 
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3.  Effect  of  Tank  Expansion 

The  flow  velocity  in  the  vent  is  very  small  in  the  beginning  compared 
to  the  98  ft/sec  characteristic  velocity  (Q/A) , expected  for  an  inviscid  fluid. 
Even  when  the  total  pressure  drop  has  increased  to  2 atm  the  mean  flow  velo- 
city is  still  only  about  25%  of  characteristic  velocity.  This  Indicates  that 
there  is  considerable  accumulation  of  liquid  in  the  tank  due  to  tank  flexing. 

4.  Effect  of  Loading  Rate 

Comparison  of  Figure  V-4  vs.  Figure  V-5  shows  that  the  rate  of  pressure 
rise  is  higher  when  the  liquid  inflow  rate  is  higher.  If  the  inflow  rate 
is  doubled,  a given  tank  pressure  is  reached  in  about  half  the  time  (time 
being  counted  after  the  tank  is  full  of  liquid). 

5.  Magnitude  of  Pressure  Rise 

From  the  results  shown  in  Figure  V-5,  it  is  seen  that  for  a two-hour 
fill  time  (which  is  an  excessive  rate  according  to  Table  lI-3)»  the  tank  pres- 
sure reaches  10  pslg  in  about  10  seconds  and  continues  to  increase.  Thus, 
tank  failure  is  expected  in  a very  short  interval  for  an  overfill  at  high 
loading  rates. 

For  a loading  rate  on  the  high  end  of  the  "normal"  range  (4-hr  time;  see 
Figure  V-4),  the  pressure  reaches  10  psig  in  about  30  seconds. 

Calculations  were  repeated  on  the  assumption  that  the  vent  is  blocked. 

The  following  equation  is  obtained  from  Equation  (21)  by  substituting 
U = 0: 


P - P„ 


(B  + <) 


in 


”i 


(26) 


Using  the  above  formula  with  appropriate  values  for  the  parameters,  it  is 
seen  that  for  a tank  being  filled  with  liquid  in  two  hours,  the  time  to 
reach  6 psig  is  between  3 and  10  seconds  which  agrees  closely  with  the  case 
of  4f  (L/D)  * 10.  Further  results  of  calculations  with  the  assumption  of 
blocked  vent  after  the  tank  is  completely  liquid  full  are  shown  in  Table  V-2. 
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Table  V-2 

Tank  Pressure  as  a Function  of  Time 
When  the  Vent  is  Blocked 
(Two  Hour  Tank  Fill  Time) 


Time 

(Seconds) 

Pressure 

(psig) 

0.5 

0.44 

1.0 

1.03 

1.5 

1.47 

2.0 

1.91 

2.5 

2.35 

3.0 

2.94 

10 

9.7 

50 

48.1 

95.84 

190.37 


The  above  comparison  indicates  that  for  these  calculations,  the  blocked  vent 
result  will  give  a reasonably  close  (and  certainly  a conservative)  estimate 
of  the  time  of  bursting  of  the  tank  shell. 


E.  Summary  of  Findings  and  Measures  to  Avoid  Overfill 


Mathematical  models  have  been  developed  to  obtain  the  rise  in  the 
tank  pressure  of  a tank  when  it  starts  venting  liquid  after  attaining  the 
liquid  full  condition.  It  is  assumed  in  these  models  that  there  is  a con- 
stant Inflow  (rate)  of  liquid  maintained.  The  models  take  into  considera- 
tion the  volumetric  expansion  of  tanks,  the  flow  friction  in  the  vent,  and 
the  pressure  drop  of  the  valves  in  the  vent  and  hydrostatic  heads  due  to 
the  vertical  sections  of  the  vent  pipe. 

The  procedure  for  determining  if  a given  loading  rate  and  vent  system 
will  present  an  overpressure  hazard  for  an  individual  tank  upon  overfill 
is  as  follows: 

2 

1.  Determine  the  venting  area  A (ft  ) from  the  cross  sectional  area 
of  the  vent  pipe(s)  loading  from  the  tank. 

. 3 

2.  Determine  the  loading  rate  Q (ft  /sec)  for  the  specific  tank  from 
the  total  pumping  rate  and  the  number  of  tanks  to  be  loaded  simul- 
taneously. 

3.  Calculate  the  quotient  Q/A  (ft/sec)  which  is  the  characteristic 
venting  velocity. 

4.  Determine  the  effective  length  to  diameter  ratio  L/D  of  the  vent 
system  using  Table  II-2  and  the  procedure  illustrated  in  Figure 
II-6. 

5.  Use  Equation  (22)  to  find  the  maximum  tank  pressure  to  be  expected 
during  cargo  overfill,  using  the  values  of  Q/A  = U and  L/D  = y/D 
obtained  in  steps  (3)  and  (4),  respectively. 

It  is  seen  from  our  results  that  when  the  loading  rate  is  excessive 
for  the  vent  capacity,  the  burst  pressure  of  the  tank  is  reached  in  a re- 
latively short  time  after  the  tank  goes  liquid  full.  Tank  damage  can  be 
avoided  by  an  automatic  shut-off  mechanism  which  detects  (a)  the  liquid 
level  in  the  tank  or  (b)  the  rapid  pressure  rise,  and  quickly  shuts  off 
the  inflow  of  liquid  to  the  tank.  Method  (b)  could  lead  to  line  rupture, 
however.  Alternatively,  overfill  tanks  and  spill  valves  can  be  provided. 
Prevention  of  liquid  overfill  condition  is  imperative  to  protect  the  tank. 

The  pressure  time  history  in  the  initial  stages  can  be  predicted  with 
reasonable  accuracy  by  assuming  the  vent  to  be  completely  blocked  off. 
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Tank  flexure  plays  an  Important  safety  role  In  relieving  overfill 
pressure.  If  the  tank  Is  Infinitely  rigid,  the  rate  of  pressure  rise 
In  the  tank  Is  extremely  rapid,  and  burst  pressure  will  be  attained 
within  a fraction  of  a second.  Our  calculations  show  that  the  delay 
times  between  filling  the  tank  and  reaching  burst  pressure  are  very 
short  (less  than  30  seconds). 

Preventative  measures  would  Include  flow  metering  devices  to  alert 

I 

i the  operator  of  excessive  loading  rates,  continuous  guaglng  systems  based 

I on  liquid  level  sensors  with  high  level  alarms,  remote  actuated  quick  clos- 

I ing  valves,  and  spill  valves. 

I 


NOMENCLATURE 


= cross  sectional  flow  area 


Units 


= velocity  of  sound  In  air  at  ambient  conditions 


ft/sec 


= Inside  diameter  of  vent  pipe 
= modulus  of  elasticity 

= fanning  friction  factor  in  the  vent  pipe 
= effective  friction  factor  2 f L/D 


acceleration  due  to  gravity 


ft/sec 


= thickness  of  bulkhead  panel 

= height  of  the  vertical  portion  of  the  vent  pipe 

= moment  of  inertia  of  panel  or  beam  section 

= ratio  of  evaporation  rate  to  the  loading  rate,  taken 
at  time  t^^^^ 

= length  of  bulkhead  panel 
= geometric  length  of  vent  pipe 
= effective  length  of  vent  pipe 
= effective  length  to  diameter  ratio 
= length  to  diameter  ratio  of  the  vent  pipe 


valve 


= equivalent  length  to  diameter  ratio  of  the  valve 
= dimensionless  mass  of  gas  in  the  tank  M/M(0) 

= rate  of  mass  efflux  through  the  vent 


Ibm/ sec 


= characteristic  mass  efflux  rate  of  air  = p Q 

a 


lb /sec 
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= rate  of  vapor  mass  addition  to  the  ullage  by  evapora-  Ibm/sec 
tion 

= mass  of  liquid  in  the  tank  at  any  time  Ibm 


= mass  of  liquid  in  the  tank  at  the  start  of  overfill  Ibm 


= rate  of  mass  inflow  (liquid  filling  rate) 


Ibm/ sec 


= pressure  Inside  the  tank  (at  the  top  of  the  tank) 
= atmospheric  pressure 


= equilibrium  vapor  pressure 
= gauge  pressure  Inside  the  tank 
= (p/p  ) “ dimensionless  pressure 

SL 

= dimensionless  hydrostatic  head  due  to  a diameter 
high  column  of  liquid  in  the  vent  pipe 

= dimensionless  gauge  pressure  inside  the  tank 

= volumetric  liquid  filling  rate 

= gas  constant  = .00232/molecular  weight  1 

= Re)Tiolds  number  for  flow  in  the  pipe  = ^ 

= stringer  spacing  or  width  of  bulkhead  panel 
= surface  area  of  cargo  liquid 


ft^/sec 

psl^-f t^-sec^/°R/lbm^ 


characteristic  tank  filling  time  (volume  of  tank 
divided  by  loading  rate) 


characteristic  vent  pipe  filling  time  (Chapter  V) 


= temperature 


mean  velocity  of  liquid  in  the  vent  pipe 


ft /sec 


mean  velocity  of  liquid  in  the  vent  pipe  at  the  start  ft/sec 
of  overflow 


k & 
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= mean  velocity  of  the  liquid  in  the  vent  pipe  if  the 
outflow  rate  is  equal  to  inflow  rate 

= tank  volume  not  occupied  by  liquid  at  any  instant 
of  time 

= total  tank  volume 

= U/U  = dimensionless  velocity  in  the  vent  pipe 
c 


= height  of  the  free  liquid  surface  (in  the  vent  pipe) 
above  the  tank  top  when  the  vent  pipe  is  being 
filled  up 

= total  length  of  the  wetted  section  of  the  vent  pipe 
when  the  vent  pipe  is  being  filled  with  liquid 

= characteristic  Mach  number 

= Mach  number  = fluid  velocity/local  sonic  velocity 
= coefficient  of  volume  expansion  of  tank  by  pressure 


■ -(-) 
v„Vdp  / 


= ratio  of  specific  heats  of  the  gas 
= turbulent  diffusion  coefficient 

= compressibility  of  liquid  = ^ ^ 


ft/sec 


.-1 

psi 


ft^/sec 


— X 4f  — ^ 
2 ^ P 


P TT^ 

— — = dimensionless  pressure  drop  in  a length 


of  vent  pipe  equal  to  its  diameter 


= kinematic  viscosity 


density  of  liquid 
density  of  gas 
density  of  pure  vapor 


2 

ft  /sec 


Ibm/ ft' 
Ibm/ft' 
Ibm/  ft' 
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Subscripts 

c 

e 

liq 

vap 

o 

T 

max 


= t/t  = dimensionless  time 
c 

= t,. ,,/t  = dimensionless  fill  time 

rill  c 


= characteristic  value 
= exit  condition 
= liquid 
= vapor 

= initial  value 

= iteration  at  time  x (Appendix  C) 
= maximum  (at  the  failure  limit) 


t 


throat 


’i 
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Appendix  A 

Alternate  Model  for  Adiabatic  Venting 

For  high  loading  rates  with  the  standpipe  vent  system  or  other  situa- 
tions where  the  gas  residence  time  in  the  vent  is  small,  the  heat  exchange 
between  the  vented  gas  and  the  ambient  may  be  negligible.  In  these  cases, 
the  appropriate  venting  assumption  is  adiabatic  rather  than  isothermal.  In 
order  to  consider  adiabatic  compressible  flow  with  friction,  the  venting 
system  may  be  considered  to  be  a converging  nozzle  with  a long  duct  attached 
as  shown  in  Figure  A^-l.  Shapiro^^^  has  developed  expressions  for  ob- 
taining the  flow  mach  numbers  and  the  mass  flow  through  such  a nozzle-duct 
system  with  friction.  These  expressions  are  utilized  to  obtain  the  mass  flow 
rate  through  the  tank  venting  system  for  a given  tank  pressure.  Using  this 
mass  flow  rate  and  equation  (11),  the  pressure-time  history  can  be  determined. 
These  steps  are  illustrated  below. 

Assuming  that  the  friction  factor  is  a constant  and  using  equations  de- 
(14) 

veloped  by  Shapiro  [Eq.  6.20  and  6.21,  page  167],  we  can  relate  the 
throat  Mach  number  (a^)  to  the  exit  Mach  number  (a^)  through  the  following 
expression: 


(A-1) 


= Y (4f  L/D)  as  the  friction  parameter. 


Where  we  have  defined  f 


Similarly,  using  equations  6.22  and  4.14-b  of  Shapiro  we  can  obtain 
an  expression  for  the  pressure  ratio  in  terms  of  and  a^: 


where  a and  a are  respectively  the  local  mach  numbers  at  the  throat  and 
t e 

at  the  exit  section.  The  mass  flow  rate  can  be  obtained  from  Equations 
(4.11)  and  (4.14-b)  of  Shapiro and  is  given  by  the  Nozzle  Mass  Flow 


j nation. 


where  p and  T are  respectively  the  absolute  tank  pressure  and  tank  gas 
temperature. 

Writing  in  Equation  (A-3)  in  dimensionless  form  to  conform  to  the 

definitions  of  tn  and  a in  the  nomenclature,  we  obtain 

c 


I 


converging  nozzle 


FIGURE  A-1  SCHEMATIC  REPRESENTATION  OF  THE  VENT  VALVE  FOR 

CALCULATION  OF  FLOW  UNDER  COMPRESSIBLE  ADIABATIC 
CONDITIONS 


A- 3 


Using  the  definition  of  (P  = p/p^  and  t = in  Eq.  (13),  and 

substituting  Eq.  (13)  and  (A-4)  into  the  mass  conservation  equation 
dM/dt  = we  obtain: 


which  is  subject  to  the  condition  f = 1 at  t = 0. 

Equations  (A-1) , (A-2) , and  (A-5)  are  three  equations  for  the  three 
unknowns  a^,  a^,  and  P.  is  determined  for  every  value  of  (P  by  solving 
simultaneously  equations  (A-1)  and  (A-2).  That  is,  for  a given  exit 
pressure  (which  is  generally  the  ambient  tank  pressure)  and  frictional 
coefficient,  and  guessing  a value  of  tank  pressure,  we  can  solve  equations 
(A-1)  and  (A-2)  for  and  a^.  This  procedure  is  described  in  Appendix  B. 
(Note  when  (P  = 1,  there  is  no  flow,  and  therefore,  and  = 0;  of 
course,  there  is  no  friction  either  under  these  circumstances.) 


As  in  the  Isothermal  case,  it  simplifies  the  form  of  the  equations 
to  rewrite  equation  (A-5)  as 


1 1 *t" 

with  P = 1 at  T = 0,  where  t = -In  (1-t);  i.e.,  t = 1-e 

The  above  equation  has  to  be  solved  numerically. 


Appendix  B 

Numerical  Solution  Procedure  for  Compressible  Adiabatic  Venting 

In  order  to  solve  the  governing  Equation  (A-1)  and  (A-2)  set  forth 
in  Appendix  A,  let 


P.  - Y-1  2 

^ ■ 2 “t 


^ “e' 

e 2 e 


(B-l) 


Then  equations  (A-1)  and  (A-2)  become 


(y-1)  ^ 6^  (1+3^) 


(B-2) 


ana 


$•  = (1+B^) 


-I  1/2 


6 (1+6  ) 

e e 


6j.  (1+B^) 


1.  e. , 


ff*  - ^ 


(1+6  J 
V't  / ‘ 


1/2 


6.  (1+6^) 

t e 


6 (1+6J 

e t 


(B-3) 


Hence,  from  (B-3)  we  get 


in  4>  - )ln  (1+6,) 


y/y-1 


in 


\ (l+6g) 

6g  (1+6^) 


(B-4) 


B-l 


Substituting  (B-4)  in  (B-2)  we  get 


Eqs.  (B-6)  and  (B-7)  are  two  equations  for  two  unknowns.  Solve 
iteratively  for  one  S in  terms  of  the  other  until  convergence  is  obtained. 


(C-5g) 


t = — — = characteristic  fill  time  for  the  vent  pipe 

c U 

c 


T = t/t  = dimensionless  time 
c 


(C-5h) 


vertical  length  of  vent  pipe 


total  length  of  vent  pipe 
M(0) tank  fill  time 


(e-5i) 


fill  M,  t characteristic  vent  fill  time 
1 c 


(C-5j) 


Using  the  definitions  in  Equations  C-5a  through  C-5j , we  write  Equation 
(C-1)  in  dimensionless  form  as 


^ = 1 
dx  (B  + k) 


[1  - v) 


X +T  — I V dx 
fill^ 


(C-6) 


with  P = 0 at  X = 0 


(C-7) 


(c-8) 


and  from  C-2 , we  have 


'■  ■ pipe  '•'d  + » I 5 


for  0 < C < H/L 


(C-9a) 


P = (L/D)p^pg  [Pd  H/L+  X £ V ] for  H/L  < C < 1 


(C-9b) 


(C-9c) 


Numerical  procedure  is  used  to  solve  for  ^ and  v as  functions  of  time. 


The  details  of  the  algorithm  are  gl  .-'-  ’ in  Table  C-1, 
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APPENDIX  D 


Evaluation  of  the  Initial  Liquid,  Velocity 

Differentiating  Equation  C-9a  with  respect  to  t , we  have  (using  Eq. 


C-8), 


/ ''  '*'> '' Irj 


At  T = 0,  we  have 


Lim  T 

T^O  J'vdT->0 

o 


Also,  from  Equation  C-6,  we  have 


1 

1 - V 

(6  + <) 

"fill+  ^ - I ^ 

^ n 

Therefore,  from  Equations  D-3  and  D-1, 


Lim  dP  1 Lim  ,,  , 

T 0 di  (3+K)Tt^. ,,  T ->■  0 ^ ^ ^pipe 


‘fill 


T 

. Lim  „ , dv  I 
+ n 2 A V — I V 

T 0 ,/*  dx  I 


dx 


x'"  0 ("d  ^ ^ 


There  are  three  cases  to  consider  for  the  above  equation: 


Case  1:  Lim  v = v = finite 

X -*■  0 ° 

Lim  ^ _ 

X 0 dx 


(B-1) 


(d-2) 


(D-3) 


finite 


Case  2: 


V = V 

o 


finite 


Lim- 
^ 0 

Lim  4^  = infinite 
T ->  0 

dv 

Lim  — 

^ 'T 

T 0 


d7  J " ■* 


finite 


Case  3: 


Lim  V = 0 

T ->■  0 


Li™  „ 4^  = infinite 

T 0 dx 

T 

Lim  dv  f , c • 

„ V I V dx  = finite 

X ^ 0 dx  J 


Case  1:  This  is  the  simplest  of  the  cases.  Substituting  v = in 
equation  D-4,  we  get  (in  view  of  D-2) 

(1  - V ) 9 

^ = (L/D)  . + Xv  ) V 

(B  + pipe  D o o 

The  above  is  a cubic  equation  in  v^.  It  can  be  shown  that  the  equation 
has  only  one  (positive)  real  solution.  This  can  be  obtained  by  using 
the  exact  method  of  solution  for  cubic  equations  (see  Handbook  of 
Mathematical  Functions  by  M.  Abramowitz  and  I.  Stegun,  page  7). 


Case  2: 


Let  C 


dx  Therefore,  v = — 

dx 


dv 

Now  Lim  = constant,  say,  c 

« clx 
X -*•  0 

with  |4^)  = V „ = finite  = v and  C = 0 at  x = 0 


Hence 


. /-  dv  dcN 

\ 't  0 


= c 


’I 


(1>5) 


I>-2 


(c/Vo) 


Therefore,  v = c/v  i.nC  + constant  of  integration 
o 


I 

i 

I 


However,  from  the  boundary  condition, 

V = V at  C = 0 
o 

This  cannot  be  satisfied  by  the  above  equation.  Hence,  Case  2 is  not 
a possible  case. 

Case  3:  By  an  analysis  similar  to  that  shown  in  Case  2,  it  can  be  shown 
that  this  is  also  a physically  Invalid  case.  Hence,  the  only  possible 
solution  to  Equation  D-4  is  given  by  the  solution  of  the  cubic  Equation 
D-5. 

Hence,  U = v x U (D-6) 

’ o o c 


I 

I 
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APPENDIX  E 

Volume  Expansion  Due  to  Deflection 
of  Bulkhead  Panels  and  Stiffeners 

Objective 

In  the  analysis  of  the  cargo  overfill  problem,  it  was  found  that  as 
the  tank  pressure  began  to  build  up  due  to  the  flow  resistance  in  the  vent 
system,  the  tank  walls  expand.  The  value  of  the  linear  expansion  parameter, 
6 = (1/V)  dV/dP,  which  is  the  fractional  volume  change  per  unit  pressure 
change,  is  very  important  to  the  analysis  and  is  estimated  as  follows: 

Approach 

Assuming  there  is  no  cargo  in  the  adjacent  tanks  to  counteract  the 
lateral  pressure  loadings  on  the  bulkhead,  we  envision  two  deflections  which 
can  be  analyzed  by  structural  formulae: 

• Deflection  of  panels;  and 


• Deflection  of  stiffners. 

This  situation  is  pictured  below  in  Figure  E-1. 


E- 


I 


First  we  must  derive  an  expression  for  the  volume  change  AV  as  a 
function  of  the  centerline  deflection  y^,  which  depends  on  the  contour 
of  the  deflection  and  the  geometry  of  the  panel.  In  addition,  we  can  use 
structural  formulas  for  y^  as  a function  of  the  pressure.  Mathematically, 
we  write  the  two  relationships  we  seek  as  follows: 


V = + AV  (y^) 


(E-1) 


y*  = f (P) 


(E-2) 


Once  we  have  these  relationships,  we  can  substitute  (E-2)  into  (E-1)  and 
differentiate  with  respect  to  pressure  to  obtain  the  expansion  parameter 


R = i ^ 

dP 


(E-3) 


Volume  Change  as  a Function  of  Centerline  Deflection 

We  take  the  panels  as  two  dimensional  since  they  are  usually  long 
rectangles.  In  general,  AV  for  an  array  of  N two-dimensional  panels 
fixed  on  both  edges  can  be  expressed  as  an  integral  over  the  panel  width 


coordinate  x: 


b 1 

AV  = Nil/  ydx  = nis  J ydn, 


(E-4) 


where  n = x/s  and  the  configuration  is  shown  in  Figure  E-1. 


k 


We  will  apply  Equation  (E-4)  to  both  elastic  and  plastic  deflection  cases. 

* 

Now  for  an  elastic  uniformly  loaded  panel,  we  have 


y = y^  16  (1  - 2n  + n^). 

elastic 


(E-5) 


This  defines  the  contour  along  the  deflected  edge.  Substituting  (E-5) 
into  (E-4)  we  obtain 


AV 


elastic 


= Ns«.  y*  16  / (n^  - 2n^  + n^)  dq 

elastic  ° 


which  gives 


'^'^elastic  = 15 


(E-6) 


elastic 


For  a plastic  uniformly  loaded  panel,  the  contour  equation  is  given 


by 


y = y*  ^ n (1  - 2n^  + n^). 
plastic 


Substituting  (E-7)  into  (E-4)  we  obtain 


16 


(E-7) 


3 . 4, 


^%lastic  - y*  . 5 / (n  - 2n  + n ) dn 

^ . o 

plastic 


A 


Bleish,  Formulas  for  Stress  and  Strain. 


which  gives 


AV  = — Ns  II  y . 

plastic  25 

plastic 


(E-8) 


Centerline  Deflection  as  a Function  of  Pressure 

For  the  plastic  and  elastic  cases  we  have  the  following  expressions 
for  a uniformly  loaded  plate: 


elastic 


= ^ _W  / ^ 

24  El  \16/ 


= -1  _W  1 \ 
* 24  El  \ 16  / 


plastic 


where  W = Load  in  pounds 

E = Modulus  of  elasticity,  and 

3 

I = Moment  of  inertia  of  panel  or  beam  section  = h £/12 
Note  that  'S ^ is  five  times  y^ 

plastic  elastic 

Now  the  pressure  P = W/sJl  so  we  have 


•ie  O /. 


1 VI 


elastic 


plastic 


24  El 


_1  M 
24  El  16 


(E-9) 


(E-10) 


Calculation  of  Expansion  Parameter  B 

Substituting  (E-9)  into  (E-6)  and  (E-10)  into  (E-8)  we  obtain 


1 NHs^P 
Aplastic  120 


(E-11) 


'rV-rr 


and 


AV 


elastic 


1 

720 


NJts^P 

El 


(E-12) 


Employing  Equation  (E-3)  we  obtain 


plastic 


N s 9} 

120  V^EI 


(E-13) 


6 


elastic 


N 

720  V^EI 


y 


(E-14) 


where  is  the  normal  tank  volume  and  N is  the  number  of  panels. 

Substituting  typical  values  for  a tank  having  four  hundred  9 x 3 ft 
panels  we  obtain 

Z = 9 ft 

= 7 X 10^  ft^ 

N = 400 
s = 3 ft 

E = 2.9  X 10^  psi 
4 

I = 9.3  in  (thickness  h = 1.0  inch) 

We  obtain  for  plastic  deflection  (E-13)  of  panels  only: 

3 = . 74  X 10  ^ (psi)  ^ 

3 

This  is  equivalent  to  a total  volume  added  of  100  ft  at  the  point  of  inci- 
dent failure  (20  psig) , which  is  approximately  1/7%  of  the  normal  tank  volume. 

A similar  analysis  for  stiffeners  was  carried  out,  and  a value  of  3 
was  obtained  approximately  equal  to  the  above  value  for  panels.  Therefore, 
the  volume  expansion  coefficient  is  estimated  to  be  approximately 


S - .0001  psi  ^ 


for  a typical  vessel.  For  any  particular  vessel,  the  value  of  3 could 
range  from  .001  to  .00001  depending  on  the  number  of  panels,  panel  dimen- 
sions, and  stiffener  spacing. 


E-5 


i- 


i{j 


